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This  program  was  directed  toward  the  understanding  of  the 


basic  colloid  thruster  generation  and  acceleration  process 
and  the  cause  of  performance  degradation.  This  thruster 
is  applicable  for  attitude  control  and  stationkeeping  on 
future  Air  Force  satellite  missions.  A theoretical  inves- 
tigation was  made  of  the  process  of  droplet  emission  to 
determine  the  source  of  negative  ion  production  causing 
sputter  erosion  of  the  emitter.  A model  of  emitter  tip 
erosion  has  been  developed  that  correlates  with  perform- 


ance degradation.  The  erosion  is  caused  by  back  bombard- 
ment of  the  emitter  by  negative  ions  produced  from  the 
breakup  of  flow  droplets  as  a result  of  collisions  with 
fast  ions.  The  calculations  were  initially  made  by  assum- 
ing that  the  droplets  were  all  the  same  size  rather  than 

using  the  distribution  of  droplet  sizes  actually  observed.  

If  at  least  one  negative  ion  is  produced  at  each  ion/  ^ 

droplet  collision,  the  resulting  sputtering  rate  is  suffi- 
cient to  explain  the  measured  erosion.  A distribution  of 
droplet  sizes  was  also  included  in  the  analysis  which 
consisted  of  a normal  distribution  of  the  charge-to-mass 
ratios.  Results  of  the  analysis  show  that  the  calculated 
collision  ratio  increases  slightly  when  the  droplet  popu- 
lation is  changed  slightly  from  uniform  size  to  the 
assumed  distribution.  The  consequence  of  varying  the 
distribution  were  examined  by  changing  the  standard  devi- 
ation and  by  replacing  the  normal  distribution  by  an 
unsymmetrical  distribution.  Recommendations  are  made  to 
eliminate  or  decrease  the  erosion  through  changing  thruster 
operational  parameters.  In  addition,  nevi;  insight  ii)to 
droplet  properties  and  collision  processes  was  obtained 
through  this  investigation. 
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Colloid  thrusters  having  multiple  needle  emitters  experience  degrada-  | 

tlon  in  performance  after  thousands  of  hours  of  operation.  This  was  i 

attributed  to  the  change  in  the  electric  field  at  the  emitter  tip 

I ' 

caused  by  erosion  which  occurred  during  that  period  of  operation. 

This  report  describes  the  development  of  an  analytical  model  for 

understanding  the  collision  processes  in  the  colloid  thruster  exhaust  j 

beam  and  emitter  erosion  mechanisms.^  Erosion  or  surface  damage  by  i 

sputtering  mechanisms  is  considered  to  be  the  cause  of  performance  ! 

degradation.  Other  erosive  (or  corrosive)  mechanisms  were  considered  | 

Initially  as  possibly  being  responsible  for  the  damage  in  the  EHD  | 

spraying  process,  but  sputtering  was  determined  to  be  the  primary  j 

i 

cause  for  producing  most  of  the  surface  damage.  j 

I 

I 

The  technical  approach  taken  initially  in  this  study  was  to  consider 

particle  interactions  and  formulate  simplified  models  involving  par-  I 

tide  collisions.  This  approach  was  used  for  the  purpose  of  gaining  i 

physical  insight  into  the  interaction  process  which  was  then  applied  I 

to  expanded  models  involving  more  complex  analyses.  ■ 

I 

The  guidelines  for  the  calculations  presented  in  this  report  were  j 

obtained  from  needle  erosion  rates  and  data.  Many  of  the  values  for 

parameters  used  in  the  theoretical  development  were  derived  from  TRW  ■' 

thruster  module  operational  data.  j 


Initially,  dlfflerent  regions  in  the  thruster  and  vacuum  chamber  were 
examined,  along  with  other  mechanisms,  to  determine  the  cause  of  the 
erosion.  The  emphasis  for  study  was  chosen  to  be  the  acceleration 


region,  near  the  droplet  emitter.  This  region  was  chosen  for  investi- 
gation for  two  reasons.  First,  the  particle  densities,  with  the  pos- 
sible exception  of  electrons,  are  highest  here  and  are  decreasing 
functions  with  distance  from  the  emitter.  Second,  negative  particles, 
suspected  as  a primary  cause  of  sputter  erosion  if  formed  as  a result 
of  collisions  in  the  acceleration  region,  are  not  generally  able  to 
escape  the  region  and  can  readily  bombard  the  emitter  surface. 

The  charged  droplet  beam  generated  by  colloid  thrusters  consists  of  a 
distribution  of  particle  sizes  covering  a range  from  molecular  ions  to 
larger  droplets  containing  as  many  as  10^  molecules/droplet.  Initially, 
this  distribution  was  supplanted  by  a typical  particle  taken  as  the 
mean  charge -to -mass  ratio,  nils  droplet  has  a radius  of  5oX  and  a 
mean  specific  charge  of  10  C/kg.  Many  aspects  of  thruster  beam 
interaction  collisions  using  this  single  particle  model  do  not  differ 
greatly  from  results  obtained  by  taking  the  actual  particle  distribu- 
tions into  account.  This  was  verified  by  extending  the  single  par- 
ticle model  to  include  charged  particle  distribution  functions. 

The  analysis,  made  in  consideration  of  the  dynamic  and  static  condi- 
tions occurring  in  the  operation  of  colloid  thruster,  supports  the 
erosion  model  presented.  The  model  involves  the  collision  of  positive 
ions  and  droplets  both  generated  at  the  emitter  and  the  subsequent 
breakup  of  the  droplets.  Among  the  fragments  produced  by  the  breakup 
of  each  droplet,  only  one  negative  ion  (probably  that  of  iodine  used 
as  the  propellant  dopant)  is  required  to  explain  the  measured  erosion. 
The  negative  ion  is  accelerated  toward  the  emitter  hitting  only 
exposed  surfaces  to  form  the  observed  erosion. 

As  a result  of  the  analysis,  recommendations  are  made  to  eliminate, 
or  at  least  minimize,  the  erosion  through  changing  thruster  opera- 
tional parameters.  A summary  of  the  erosion  model  resulting  from  this 
program  is  given  in  subsection  7.1,  xjhich  is  recommended  for  review. 
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SECTION  2 


THKUSTER  EMITTER  DYNAMICS 


2.1  PERFORMANCE  DEGRADATION  AND  EROSION 


Colloid  charged  particle  generation  was  easily  accomplished  in  the 
initial  stages  of  the  research,  but  the  control  and  understanding  of 
the  fundamental  physical  processes  involved  proved  a more  difficult 
task.  The  development  of  colloid  thrusters  has  proceeded  primarily 
by  empirical  Investigations  without  sufficient  and  practical  theor- 
etical models  to  describe  the  interdependence  of  such  thruster 
variables  as  thrust,  specific  impulse,  and  specific  charge  of  efficiency. 

Duration  tests  of  colloid  thruster  modules  at  TRW,  Inc.,  show  degrada- 
tion in  performance  after  only  several  thousand  hours  of  operation. 

Examination  of  the  emitters  revealed  some  form  of  erosive  attack  at 

2 

the  capillary  needle  edges.  The  performance  degradation  measured  in 
terms  of  decreasing  specific  Impulse  has  been  attributed  to  this  dam- 
age. The  damage  is  thought  to  be  material  sputtering  erosion  caused 
by  negative  particle  bombardment.  However,  because  of  the  extraordi- 
nary conditions  at  the  emitter  tips,  other  surface  damage  mechanisms 
may  be  active  which  would  not  be  anticipated  in  a less  extreme  environ- 
ment. Data  from  long  term  operation  of  emitters  reveal  erosion  that 
occurs  at  the  rate  of  about  one  mil  in  a year.  The  emitter  tip 
geometry  is  apparently  altered  sufficiently  to  modify  the  electric 
field  in  the  emission  areas  and  thereby  decrease  the  specific  impulse. 
Performance  can  be  restored  by  increasing  the  high  voltage  applied  to 
the  emitters,  but  this  does  not  completely  restore  the  local  electric 
field  to  reproduce  the  desired  operation.  A somewiiat  different  oper- 
ating mode  is  thereby  achieved  in  which  degradation  may  be  accelerated, 
in  addition  to  other  unwanted  effects  such  as  changes  in  the  beam 
spread  and  beam  distributions. 
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A photograph  made  by  TRW,  Inc,,  showing  the  typical  pattern  of  emitter 
damage  is  shown  in  Figure  1 (Ref.  2),  The  individual  sites  for  jet 
formation  can  be  observed,  along  with  marks  due  to  machining  made 
during  the  emitter  fabrication.  The  rims  are  eroded  only  on  the  out- 
side, indicating  that  attack  by  electrolytic  action  is  not  prevelant. 
The  physical  evidence  of  erosion  on  the  rims  and  outside  needle  diam- 
eter suggests  the  damage  is  caused  by  negative  ion  sputtering  with 
the  scalloped  mounds  protected  by  the  liquid  jets.  Figure  2 illus- 
trates the  eroded  area  of  the  needle,  showing  significant  dimensions. 
The  results  from  one -dimensional  modeling  of  the  interaction  phenomenon 
in  the  colloid  exhaust  plume  is  advanced  as  an  explanation  for  the 
erosive  damage  to  colloid  thruster  emitters. 

2.2  THE  ACCELERATION  REGION 


The  most  important  area  of  interest  in  the  erosion  study  is  the  accel- 
eration region  in  the  vicinity  of  the  emitter  tip.  Figure  3 depicts 
the  TRW  needle  geometry  showing  the  acceleration  region,  or  region  of 
interaction,  under  consideration  for  particle  collisions  in  the  sput- 
tering studies.  The  needle  emitter  has  a 0.014  inch  outer  diameter 
and  a 0.005  inch  inner  diameter  and  is  constructed  from  a platinum 
(907o)/irtdium  (107.)  tubing.  The  needle  is  surrounded  by  a concentric 
cylindrical  shield  having  a 0.062  inch  outer  diameter.  The  shield/ 
emitter  assembly  is  centered  in  a planar  extractor  electrode  which 
has  a 1/8  inch  aperture. 

The  acceleration  region  in  Figure  3 extends  from  the  needle  emitter 
tip  to  a distance  approximately  equal  to  the  radius  of  the  extractor 
aperture.  It  includes  the  particle  beam,  diverging  from  the  emitter, 
and  the  intense  electric  field. 


The  Potential  Function 


An  analytical  expression  for  the  potential  distribution  arising  from 
the  electric  fields  in  the  region  of  the  emitter  and  extractor  elec- 
trodes was  formulated.  The  potential  function  is  required  for  calcu- 
lations of  charged  particle  density  functions,  cross  sections,  and 
other  analyses  involving  charged  particles.  The  required  distribution 
function  would  be  ideally  derived  from  the  boundary  conditions  implied 
in  electrode  geometry  shown-  in  Figure  3 and  also  include  the  effects 
of  space  charge.  A potentjial  function  found  to  be  adequate  for  much 
of  the  quantitative  analyses  in  this  report  is  given  below; 


In 


V(x)  = 


,x  + r r 
2 o> 


f 


(1) 


This  expression  was  deijived  by  considering  a tip  geometry  of  parabolic 
shape  surrounded  by  a grounded  surface  having  a parabolic  shape,  where 
V(x)  and  are  the  potentials  at  a distance,  x,  from  the  emitter  and 
at  the  emitter  surface,  (respectively.  The  distance  between  the  tip 
and  the  grounded  surface' given  by  R and  r is  the  radius  of  curvature 
at  the  tip.  For  application  of  Eq.  (1)  to  the  emitter  geometry,  R is 
taken  as  the  radius  of  the  extractor  electrode  aperture  and  r^  as  the 
radius  of  the  needle  emitter.  For  the  case  where  a needle  emitter  is 

i 

surrounded  by  a shield  el^trode,^  shown  in, Figure  3,  r^  is  taken  as 
the  radius  of  the  shield  electr^e.  The  pbtential  calculated  using 
Eq.  (1)  falls  to  a value  of  zeijb  atl  a distance,  R,  along  the  emitter 
axis,  equivalent  to  the  extraetpr  aperture  radius.  The  value  of  r 
determines  the  rate  at  which  the\  field  deert^ases  with  distance. 

Values  for  the  potential  distributiort  «(^re  calculated  for  the  geometry 
used  for  computer  generated  field  plots  ^iipplied  by  TRW,  Inc.^  The 
results  are  shown  :^n  Figure  4,  where  curve  A was  calculated  from  the 


f 


potential  function  (Eq.  1)  and  curve  C was  obtained  from  the  computer 

studies.  The  differences  between  the  computer  data  and  calculations 

can  be  seen  from  examination  of  Figure  4.  At  x = R - 1/2  r , the 

o 

potential  functions  goes  to  zero,  whereas  the  computer  values  extend 
slightly  beyond  R.  Compared  with  the  analytic  expression,  the  com- 
puter potential  shows  a more  rapid  decrease  as  a function  of  distance 
from  the  emitter  tip.  Curve  B is  a plot  of  Eq.  (1)  f^r  the  geometry 
where  erosion  data  are  not  available.  The  needle  voltage  is  assumed 
to  be  13  kV,  and  r^  is  the  radius  of  the  shield  equal  to  0.031  inch. 
Note  in  Curve  B that  the  field  does  not  drop  off  as  rapidly  in  the 
vicintiy  of  the  emitter  tip  as  a result  of  the  presence  of  the  shield 
electrode. 

2.3  SPUTTERING  AS  A MECHANISM 

Emitter  surface  sputtering  rates  are  examined  which  can  account  for 
the  loss  of  material  observed  experimentally  during  long-term  operation 
of  colloid  thrusters.  For  every  ion  that  bombards  a surface,  S atoms 
are  sputtered  off  the  surface,  where  S is  no  more  than  an  order  of 
magnitude  from  unity  for  the  species  and  energies  involved  here. 

The  bombarding  ion  current,  I,  necessary  to  sputter  a mass,  m,  of 
target  material  is  given  by  the  relation: 

1.66  X lo’-^  V p 

I (ions/sec)  = ^ g , 

where,  V is  the  sputtered  volume,  p the  density  of  the  material 
s 

(pg/unit  volume),  A is  the  mass  number  of  the  target  atom,  S is  the 
sputtering  yield  (a  oms/ion),  and  t is  the  bombardment  time  in  hours. 

In  order  to  calculate  the  incident  ion  bombardment  rate  at  colloid 
emitters  that  could  cause  erosion,  assuming  various  sputtering  yields. 
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published  erosion  rates  for  the  needle  emitters  tested  at  TRW  were 
used.^  Assuming  a 0.001  inch  thickness  for  the  needle  rims  and  an 
erosion  rate  of  1 mil/year,  the  following  incident  negative  ion  rates 
were  calculated: 


Sputtering  Yield 
S (atoms/ion) 

1 

5 

10 


Incident  Ion  Rate 
I (ions/sec) 

1.39  X 10^ 

2.78  X 10® 

1.39  X 10® 
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Bombardment  rates  on  the  order  of  10  to  10  ions/sec  could  account 
for  the  observed  erosion  rates.  Sputtering  yields  in  the  range  from 
1 to  10  atoms/ion  are  possible  if  the  bombarding  species  is  assumed  to 
be  negative  iodine  ions  or  large  negative  cluster  ions.  To  our  knowl- 
edge, there  are  no  data  available  on  the  sputtering  yields  for  negative 
iodine  atoms  incident  on  a platinum  target.  A crude  approximation  for 
the  sputtering  yields  involving  iodine  ions  can  be  made  by  extrapola- 
ting data  and  theory  for  other  species  incident  on  platinum.  Figure  5 
shows  a plot  of  sputtering  rates  for  argon  ions  on  platinum  over  the 
energy  range  1 to  13  kV.  The  argon  yield  curve,  denoted  by  S(A^)  was 
constructed  using  a universal  yield-energy  relation  for  argon  ion 
bombardment  on  various  solid  targets.  The  upper  curve  in  Figure  5 is 
the  estimated  sputtering  ratio  for  iodine  ions  incident  on  platinum 
denoted  by  S (A+)  (Mj^/M2).  This  curve  was  constructed  from  the  argon 
curve  assuming  that  sputtering  yields  on  a given  target  by  various  ion 
species  is  directly  proportional  to  the  mass  of  the  incident  ion. 

This  was  obtained  by  multiplying  the  sputtering  ratios  for  argon  by 
the  ratio  of  the  iodine  to  argon  ion  masses  which  is  3.14.  The  result- 
ant yield  curve  for  iodine  is  believed  to  overestimate  the  sputtering 
ratio,  and  the  more  realistic  sputtering  ratios  probably  fall  between 
the  argon  and  the  estimated  iodine  yield  curves.  Therefore,  it  is  not 
unreasonable  to  assume  that  iodine  ions  could  sputter  as  many  as 
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SPUTTERING  RATIO  (atoms/ion) 


TaSET! 


10  surface  platinum  atoms  percollision  for  the  energy  range  considered. 
If  negative  Ions  of  the  type  OH"  or  negatively  charged  glycerol  frag- 
ment Ions  are  also  contributing  to  the  erosion  of  emitters,  they  are 
not  expected  to  produce  high  sputtering  yields  similar  to  iodine,  since 
their  masses  are  much  smaller. 
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SECTION  3 


PARTICLE  SPECIES  IN  EXHAUST  BEAM  ACCELERATION  REGION 


Before  collision  rates  and  cross  section  analyses  can  be  undertaken, 
determinations  will  be  made  on  particle  species  populating  the  exhaust 
beam  and  their  corresponding  density  functions.  The  particles  most 
likely  involved  in  beam  collision  processes  are  thermally  evaporated 
neutral  glycerol  molecules,  EHD  generated  molecular  ions  and  charged 
droplets,  and,  finally,  electrons.  With  the  exception  of  the  latter, 
approximate  density  functions  can  be  determined  with  some  knowledge 
of  the  generating  mechanisms  producing  the  particles.  Since,  in  this 
analysis,  we  are  not  concerned  with  particle  densities  outside  the 
accelerating  region,  variations  in  charged  specie  densities  will  be 
due  to  particle  accelerations  and  beam  divergence.  The  density  func- 
tion for  neutral  glycerol  molecules  was  calculated  starting  with  a 
conventional  cosine  distribution. 

3.1  THERMALLY  EVAPORATED  NEUTRAL  GLYCEROL  MOLECULES 

The  number  of  molecules  per  square  centimeter  per  second  (J)  passing  a 
unit  area  per  unit  time  at  a distance,  r,  from  a source  whose  area 


is  A is: 


O V A cosG 
*o  o 

- 2 
4 n r 


where  p is  the  number  of  molecules  per  cubic  centimeter  at  the  source 
_ o 7 

and  V is  the  mean  velocity  of  the  molecules.  Equation  (2)  applies 
o 

when  r is  large  compared  with  the  dimensions  of  the  source.  In  the 
accelerating  region,  distances  are  not  large  compared  to  the  needle 
dimensions.  A point  source  approximation  such  as  Eq.  (2)  is  not  valid 
for  accurate  determinations  of  the  density  of  thermal  glycerol  mole- 
cules in  close  proximity  to  the  evaporating  surface  of  the  propellant 
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fluid.  The  problem  of  an  extended  effusive  source  was  solved  for  the 
special  case  of  determining  the  number  of  density  along  the  colloid 
emitter  axis  perpendicular  to  the  emitter  surface  (0  = 0°).  The  result 
is  given  by: 


J(x)  = 


where  x is  the  distance  along  the  emitter  axis  and  r^  is  the  radius  of 

the  emitter.  At  x = 0,  the  expression  for  J(x)  reduces  to  k:  v^, 

where  p is  the  density  at  the  surface  and  v is  the  mean  thermal  veloc- 
o o 

ity.  When  x is  large  compared  with  r^  and  for  the  case  0 = 0°,  Eq.  (3) 
reduces  to 

J(x)  = 


which  is  identical  to  Eq.  (2)  for  a source  having  circular  dimensions 
2 

of  nr  , A normalized  form  of  Eq.  (3),  i.e.,  f(x)  is  graphically  dis- 
o 

played  as  a function  of  distance  from  the  source,  x,  in  units  of  x/r 

o 

in  Figure  6,  For  comparative  purposes,  the  point  source  approximation 
given  by  Eq.  (4)  is  also  plotted.  At  large  distances  from  the  emitter, 
the  extended  source  solution  can  be  seen  to  approach  the  point  source 
formulation.  When  x is  small,  the  point  source  solution  diverges  to 
infinity  while  the  normalized  flux  approaches  unity.  Note  that,  at  a 
distance  of  x = 10  r^,  the  particle  density  will  have  decreased  by  a 
factor  of  100. 

The  neutral  glycerol  molecular  density  as  a function  of  x is  defined  as: 

p(x)  = J(x)/7 
o 
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Variations  In  the  density  of  thermal  glycerol  molecules  in  the  accel- 
erating region  can  be  calculated  using  Eq.  (8),  remembering  that 

and  V are  temperature  dependent.  At  a temperature  of  20°C,  the  vapor 
^ "A 

pressure  of  glycerol  is  1.7  x 10  torr.  The  calculated  value  of  J 

16  2 ^ 
from  Eq.  (6)  is  given  by  3.7  x 10  mol/cm  -sec.  Using  this  value 

and  the  velocity  given  by  Eq.  (7),  the  neutral  glycerine  molecular 

density  given  by  Eq . (8),  at  a distance  of  8.75  mils,  is  given  by 
11  3 

2 X 10  mol/cm  . At  10  emitter  radii,  the  density  drops  to  about 
9 3 

2 X 10  mol/cm  and  decreases  rapidly  with  greater  distance. 

3.2  MOLECULAR  ION  CURRENT  DENSITIES 

Molecular  ions  are  defined  as  any  cluster  composed  of  glycerol  mole- 
cules that  carry  a single  electronic  charge.  The  existence  of  molecular 
ions  is  evidenced  by  the  peaks  occurring  in  the  Initial  portion  of  time- 
of -flight  waveforms,  and  they  have  also  been  detected  using  EH  ioniza- 
tion mass  spectrometry  techniques  employed  to  measure  the  ionization  of 
glycerol  from  needle  emitters. In  the  latter  study,  singly  charged 
ions  containing  up  to  eight  glycerol  molecules  were  observed.  The 
charge  associated  with  the  molecules  of  the  form  [g  + Na]^  for  n from 
0 to  8 was  attributed  to  the  Na  ion.  During  typical  operation  of  a 
colloid  needle  emitter,  10  percent  of  the  total  needle  current  can  be 
attributed  to  ion  emission. 

Data  taken  from  Reference  10  was  combined  with  the  ion  peak  normally 
observed  in  time-of -flight  (TOP)  curves  to  produce  Figure  7.  The 
dashed  line  is  a typical  ion  peak  in  TOP  data  with  the  total  ion  cur- 
rent plotted  on  an  arbitrary  scale  of  100.  The  different  molecular 
ion  fragments  taken  from  Reference  10  were  plotted  as  vertical  lines 
at  the  time-of-f light  for  the  arrival  of  such  ions  as  Na^,  G + Na^ 

(where  G represents  a glycerine  molecule),  etc.  Since  the  TOP  data 
have  a relatively  long  time  constant,  the  scope  trace  was  not  able  to 
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Figure  7.  TOF  Construction  of  Ion  Peak  Using  EHD  Mass  Spectrometry  Data 


follow  the  sharp  steps  shown  with  the  continuous  line  but  showed  the 
decay  seen  in  the  daslied  line.  Tlie  correlation  of  these  two  data 
verify  the  existence  of  these  various  ion  fragments  in  the  colloid  beam. 

The  emission  rate  of  ions  can  be  calculated  from 


where  I.  is  that  portion  of  the  needle  current  attributed  to  ion  emis- 

-19 

Sion  and  q is  the  charge  of  the  ion  equivalent  to  1.6  x 10  coulomb. 
Assuming  a total  needle  current  of  10  pA,  then  1 pA,  and  molec- 

ular ion  emission  rates  are  calculated  to  be  approximately  6 x 10^^  ions/ 
sec.  It  is  interesting  to  note  that  ion  emission  -ates  are  of  the  same 
order  of  magnitude  as  droplet  emission  rates  to  be  discussed  next. 


3.3  CHARGED  DROPLETS 

A number  density  for  a charged  droplet  of  given  size  in  the  accelerating 
region  can  bo  calculated  from  the  particle  generation  rate  at  the  emit- 
ter. The  charged  particle  flux  0,  in  droplets  per  second  per  unit  area, 

is  determined  from  0 - r/A  where  A is  the  emission  area  at  the  surface 

r r 

of  the  omitti^r  and  U is  the  emission  rate  in  droplets  per  second.  Assum- 
ing the  colloid  beam  consists  of  particles  having  uniform  size  determined 
from  the  mean  charge-to-mass  ratio  of  the  distribution,  a characteristic 
frequency  for  this  particle  over  the  entire  emitter  is  defined  as 

V = I^/q,  (9) 

where  is  the  total  needle  current  and  q is  the  charge  carried  by  the 
mean  specific  cliargi  particle.  Applying  Rayleigh's  criteria  determines 
the  charge  (q)  for  the  moan  specific  charge  particle  for  glycerol  and 
the  characteristic  emission  frequency  becomes 

u = 1.49  X 10^^  <C>  , 
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where  (C)  is  the  mean  specific  cliarge  characterizing  the  distribution. 

' 4 , 12 

For  a mean  specific  charge  of  10  C/kg  and  of  10  \i.\,  v = 1.49  x 10 

droplets/second . The  mean  specific  charge  is  given  by  (c)  = I^/m,  where 

ra  is  generally  one  p,g/sec  and  I is  given  above.  This  results  in  the 

A 

value  of  10  C/kg. 


The  charged  particle  flux  density  in  the  accelerating  region  is  given  by 
the  continuity  relationship  as 

i)  - ^(x)  v(x)  (10) 

The  velocity  ’'(x)  of  the  droplet  is  not  constant  because  of  the  accel- 
erating forces  and  varies  according  to 

v(x)  = 2C[V^  - V(x)  ] (11) 

wlierc  and  V(x)  are  the  potentials  at  tlie  emitter  surface  and  at  a 
distance  x from  the  cmitLei\,  respectively.  Assuming  p(x)  v(x)  to  be 
constant  neglects  any  divergence  due  to  off-axis  particle  emission 
and  space  cliarge  spreading.  Combining  Eqs . (10)  and  (11),  the  droplet 
density  function  becomes 

P(x)  = 0 1^2  C [V^  - V(x)1^J  (12) 

Note  that  V(x)  is  the  potential  functijjn  in  Eq.  (1)  of  Section  2.  This 
solution  for  the  density  function  will  diverge  in  the  limiting  case  at 
X = 0 where  the  velocity  approaches  zero  and  p(x)  -•  This  situation 
is  commonly  encountered  in  problems  of  space  charge  limited  currents. 


The  droplet  density  function  as  a function  of  distance  from  the  emitter 
taken  from  Eq . (12)  is  shown  in  figures.  For  this  plot  (c)  was  taken 
as  10^  C/kg,  as  10  uA,  as  13  kV,  and  the  standard  geometry  shown 

in  figure  3 was  used  > calculate  values  for  the  potential  function.  A 


seen  in  the  figure,  droplet  densities  are  on  the  order  of  10  per  cubic 
centimeter  in  the  accelerating  region.  Because  of  the  space  charge 
divergence,  the  density  will  decrease  more  rapidly  than  indicated  here. 

3.4  STRUCTURE  AND  ENERGETICS  OF  CHARGED  DROPLETS 

The  charged  droplet  distribution  represented  by  the  mean  specific  charge 

4 

of  10  C/kg  has  a significant  role  in  the  collision  studies.  The  pre- 
diction of  collision  products  after  a droplet  is  bombarded  by  another 
particle  requires  some  knowledge  of  the  structure  and  energy  of  the 
droplet.  An  examination  of  particle  properties  including  a discussion 
of  particle  stability  and  energetics  will  be  made  and  subsequently 
applied  to  collision  analyses  in  the  next  section. 

3.4.1  DROPLET  PROPERTIES 

Table  I lists  some  typical  properties  of  a droplet  assumed  spherical 
in  shape.  The  oroplet  mass,  charge,  and  radius  were  deternnned  from 
Rayleigh's  criteria  for  a specific  charge  of  10^  C/kg.  The  ratio  of 
glycerol  molecules  to  sodium  iodide  molecules  was  calculated  for  a 
solution  of  20  grams  Nal  dissolved  in  100  ml  of  glycerol.  The  ratio 
was  then  used  to  determine  the  number  of  glycerol  and  Nal  molecules  in 
the  droplet  from  the  following  relation.ships; 


N 


% 


(13) 


and 


% 


Nal 


(14) 


where  m^  and  ^re  tlie  masses  of  the  glycerol  and  Nal  molecules  in 
grams,  is  ttic  mass  of  Llie  droplet,  and  R is  the  ratio  of  glycerol 
to  Nal  molecules. 
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TABLE  1 


DROPLET  PROPERTIES 


Charge- to-Mass  Raiio 

Radius 

Charge 

Mass 

Number  of  Glycerol  Molecules 
Number  of  Nal  Molecules 
Number  of  Surface  Molecules 
Glycerol/Nal  Molecule  Ratio 
Emission  Frequency 
Surface  Area 
Volume 

Geometric  Cross  Section 


10^  C/kg 
o 
SOA 

+40e 

6.67  X 10  kg 

3.74  X 10^ 

3.7  X 10^ 

10^ 

10: 1 

10^*^  - 10^^  drops/sec 

-12  2 
3.14  X 10  cm 

-19  3 

5.24  X 10  cm 

-13  2 

7.85  X 10  cm 


Although  dissolved  iodine  atoms  contribute  little  to  the  mass  of  the 
droplet,  their  presence  in  the  droplet  may  be  a key  factor  in  the  pro- 
duction of  negative  ions  during  droplet  breakup  or  evaporative  colli- 
sions (see  Section  6).  The  droplet  emission  frequency  of  10  to  lo 
drops/sec  was  estimated  from  the  characteristic  frequency  derived  pre- 
viously. The  frequency  in  Table  I represents  the  rate  of  droplet  emis- 
sion from  a single  spraying  site  and  Implies  a number  of  individual 
spraying  sites  between  10  and  100.  Evidence  that  30  to  40  spraying 
sites  are  active  during  operation  of  needle  emitters  was  seen  in  scan- 
ning electron  photomicrographs,  Figure  1.^  A frequency  spectrum  con- 
structed from  time-of-f light  waveform  supplied  by  TRW  shows  that  a 

4 

particle  of  mean  specific  charge  of  10  C/kg  has  an  emission  frequency 
(Al/q)  in  the  range  10^^  to  10^^  drops/second. 

The  stability  of  a charged  droplet  generated  in  the  EHD  spraying  process 
will  depend  to  a large  extent  on  the  intermolecular  forces  of  attraction 
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acting  between  individual  molecules  in  tbe  droplet.  These  forces,  as 
shall  be  seen  later,  which  arc  de.scribcd  by  Van  dcr  Waals  dipole  and 
dipole- induced  interactions,  arc  not  particularly  strong  compared  with 
intramolecular  covalent  or  ionic  bonds.  Further,  the  cohesive  sta- 
bility of  a charged  droplet  is  lessened  by  the  electrostatic  repulsion 

energy  due  to  excess  charge  and  also  by  elevated'  temperatures  arising 

2 

from  1 R heating  during  the  droplet  formation  process. 

3.4.2  ENERGETICS  OF  AN  UNCHARGED  DROPLET 


Insight  into  the  stabibility  of  a charged  liquid  droplet  can  be  gained 
by  considering  the  various  energy  states  of  a molecule  residing  in  an 
uncharged  droplet.  The  energy  level  diagram  for  a liquid  droplet  is 
shown  in  Figure  9,  and  the  following  terms  arc  defined  in  reference 
to  the  figure; 


E^  = absolute  energy  of  molecule  on  surface  of  droplet 

E^  = absolute  energy  of  molecule  in  the  interior  of  droplet 
(cohesive  energy)  in  normal,  bound  state 

E = energy  state  of  unbound,  free  molecule 
r 

c = surface  energy  - energy  required  to  bring  molecule  from 
interior  of  droplet  to  the  surface 

E^  = total  energy  of  all  molecules  in  droplet  occupying 
surface  and  interior  states 


The  total  cohesive  energy  of  all  molecules  in  the  droplet  is  the  sum 
of  the  energies  of  the  molecules  in  the  droplet  interior  and  the  mole- 
cules on  the  surface: 


= (N, 


- N ) 
s 


N E 
s s 


(15) 


where  N^  and  N are  the  total  number  of  molecules  in  the  droplet  and 
t s 

the  number  of  molecules  on  the  surface,  respectively.  Rearrangement 
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of  the  terms  in  Eq . (15)  yields: 


E = N E - N (E  - K ) (16) 

ttLs'Ls*^ 

VA 

From  the  cneixy  diogrom,  E - E - e and  e can  be  replaced  by 

& s s ^ • 

where  Y is  the  surface  energy  and  A is  the  surface  area  of  the 
droplet.  Making  the  proper  subsl i Lutions,  the  total  energy  becomes; 


E = N E - yA 
t t L 


Thus,  E^  represents  the  cohesive  energy  of  a droplet  whose  value  can 
be  estimated  from  the  actual  heat  of  vaporization  when  the  substance 
evaporates  into  the  atomic  or  molecular  constituents.  The  energy,  E^, 
can  also  be  interpreted  as  the  energy  required  to  evaporate  the  droplet 
into  a gas  composed  of  its  individual  molecules  assuming  no  energy 
losses  in  the  form  of  kinetic  energy,  ionization,  excitation,  etc. 

For  the  case  of  a charged  droplet,  the  energy  will  be  less  by  an 
amount  corresponding  to  the  electrostatic  repulsive  energy  due  to 
the  mutual  interaction  of  the  net  charges  residing  on  the  surface  of 
the  droplet.  The  electrostatic  energy  of  a system  of  cliarges  on  the 
droplet  surface  whose  radius  is  R and  tlie  number  of  charges  is  Z can 
be  expressed  as: 


An  € R 
o 


The  total  cohesive  energy  of  a charged  droplet  becomes 


E„  = N E,  - yA  > f — ^ . 

t t L 4n  e R 

o 
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The  following  subsection  is  devoted  to  a discussion  of  the  energy 
of  attraction  between  molecules  in  a liquid  droplet. 

3.4.3  INTERMOLECUL^R  ATTRACTION  ENERGIES 

To  simplify  the  development  of  a model  for  droplet  energy,  the  droplet 
is  treated  as  though  composed  of  glycerol  molecules  only.  The  conse- 
quences of  neglecting  the  effects  of  Nal  molecules,  I atoms,  and  Na^, 
l"  atomic  ions  on  the  total  intermolecular  binding  energy  of  a droplet 

will  be  discussed  later.  Glycerol  molecules  possess  a permanent  dipole 

“ 18  11 
moment  of  2.68  x 10  e.s.u.  which  is  greater  than  that  of  water. 

The  forces  of  attraction  acting  between  glycerol  molecules  are  due  in 

part  to  its  dipole  moment.  The  Van  der  Waal's  intermolecular  attraction 

energy  between  two  identical  molecules  possessing  permanent  dipole 

moments  is  given  by  Reference  12. 


* 


3 2 \ 

4 ^ y ’ 


where  the  interaction  distance  between  the  molecules  is  given  by  r, 

(i  is  the  dipole  moment,  Q:  is  the  polarizability  of  the  glycerol  mole- 
cule, k is  Boltzmann's  constant,  T the  absolute  temperature,  and  I is 
the  ionization  potential  of  a glycerol  molecule.  The  first  term  inside 
the  parentheses  is  the  average  energy  of  the  dipole-dipole  interaction. 
The  second  term  represents  the  energy  of  attraction  due  to  the  induced 
dipole  in  the  'first  molecule  by  the  field  of  the  permanent  dipole  of 
the  Second  molecule,  and  vice-versa.  The  last  term  is  the  dispersion 
energy  representing  a net  force  due  to  instantaneous  forces  of  attrac- 
tion resulting  from  the  appearance  of  momentary  dipoles. 


Before  E , the  intermolecular  attraction  energy  between  glycerol 
molecules,  can  be  calculated,  the  values  of  the  polarizability  (o), 
ionization  potential  (I),  and  interaction  distance  (r)  for  glycerol 


27 
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must  be  estimated.  The  polarizability  of  the  glycerol  molecule  was 
estimated  using  the  relation: 


where  K is  the  dielectric  constant,  M is  the  molecular  weight,  p is 
the  density  of  the  liquid,  and  N is  Avogadro's  number. 


Using  Eq.  (20),  the  polarizability  for  glycerol  was  calculated  to  be 
2.72  X 10  cm^ , which  is  the  value  found  in  most  organic  liquids. 

In  general,  the  polarizability  tends  to  increase  with  increasing  size 
of  the  molecule.  Thus,  in  calculating  the  polarizability  of  glycerol, 
the  electronic  and  atomic  contributions  were  neglected  because  they  are 

UK 

small  by  comparison. 


The  polarizability  obtained  here  is  therefore  the  orientation  polariza- 
bility which  expresses  the  tendency  for  a molecule  to  align  itself  in 
an  electric  field.  It  was  mentioned  previously  that  the  ionization 
potential  of  organic  molecules  did  not  vary  over  wide  limits,  and 
inspection  of  data  for  measured  ionization  potentials  shows  that  an 
average  energy  of  10  eV  can  be  used  for  glycerol.  The  greatest  un- 
certainty involves  an  estimation  of  the  interaction  distance,  r, 
between  glycerol  molecules.  For  the  present  analyses,  r is  taken  to 
be  the  diameter  of  the  glycerol  molecule  (D)  calculated  using:^^ 

-8 

D = 1.33  X 10  { -^  ) cm  . (21) 


where  M is  the  molecular  weight  and  r is  the  density  of  the  liquid  in 
3 

g/cm  . The  results  for  the  evaluation  of  EL  are  particularly  sensitive 

o ^ 

to  the  value  of  D (5.56A)  because  of  the  r dependency  in  Eq.  (19). 


The  individual  contributions  of  the  dipole,  induced  dipole,  and  dis- 
persion energies  to  the  total  interaction  attractive  energy  E^^  in 
Eq.  (19)  as  calculated  for  glycerol  are  shown  in  Table  II. 
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TABLE  II 


CONTRIBUTIONS  (7,)  TO  VAN  DER  WAAL'S  ENERGY  FOR  GLYCEROL  MOLECULES 


Dipole 


Induced  Dipole 


Dispersion 


Thus,  we  can  see  that  the  dipole  coupled  forces  of  attraction  account 
for  about  12  percent  of  the  intermolecular  binding  energy  between  the 
molecules  in  liquid  glycerol.  The  total  Van  der  Waal's  energy  between 
two  glycerol  molecules  was  calculated  to  be  E^  = 0.2  eV.  This  Van  der 
Waal's  bonding  energy  is  on  the  order  of  10  times  the  intermolecular 

bonding  energy  between  H„  molecules  and,  in  general,  these  attractive 

^ 16 
energies  are  stronger  with  increasing  size  of  the  molecules.  The 

energies  involved  in  intermolecular  forces  between  molecules  can  also 
be  estimated  from  the  heat  of  vaporization  of  glycerol.  The  heat  of 
vaporization  of  glycerol  is  approximately  9.5  x 10^  joules/kg.  If  we 
multiply  this  value  of  by  the  mass  of  a glycerol  molecule  and  con- 
vert to  electron  volts,  we  find  the  energy  required  to  break  the  bonds 
existing  between  glycerol  molecules  is  ~ 0.9  eV.  This  value  appears 
to  be  high,  but  this  is  consistent  with  the  fact  that  glycerol  has  an 
unusually  low  vapor  pressure.  Theinterraolecular  bonding  energies  in 
glycerol  determined  from  the  heat  of  vaporization  are  4.5  times  that 
calculated  by  Van  der  Waal's  energies.  The  discrepancy  can  be  explained 
in  part  by  the  choice  for  the  interaction  distance  in  Eq.  (19),  which 
depends  on  the  assumed  diameter  of  the  molecule.  A drecrease  in  r of 
about  20  to  25  percent  in  the  Van  der  Waal's  equation  will  increase 
the  binding  energy  a factor  of  four, in  closer  agreement  with  the  heat 
of  vaporization.  Perhaps  even  more  important,  Eq.  (19)  is  based  upon 
the  binding  energy  of  two  molecules  rather  than  the  multimolecular 
forces  that  exist  in  a solid  or  liquid.  Thus,  the  number  of  bonds 
might  be  increased  m nyfold  to  account  for  this,  thereby  resulting  in 
a much  larger  value  of  Ej^.  In  view  of  the  uncertainty  in  the  calcula- 
tion, we  shall  use  the  value  of  E^^  = 0.9  eV  for  the  total  cohesive 
energy  of  the  charged  droplet. 
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The  total  energy,  E^,  of  the  molecular  droplet  system  for  a 50A  droplet 
at  20  C was  calculated  using  Eq.  (18).  Tlie  relative  contribution  for 
each  energy  terms  is  shown  below: 
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In  the  calculation,  was  taken  as  the  sum  of  the  glycerol  and  sodium 
iodide  molecules  in  the  droplet,  assuming  that  the  intermolecular  at- 
tractive forces  acting  between  two  glycerol  molecules  is  essentially 

the  same  as  the  forces  acting  between  a Nal  molecule  and  a glycerol  , ] 

S 

molecule.  Therefore,  the  total  energy  of  the  droplet  is  'V  3000  eV.  | 

In  the  section  on  particle  collisions,  the  droplet  binding  energy  i 

plays  a key  role  in  the  examination  of  probable  mechanisms  for  negative 

Ion  formation.  ‘ 


I 
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SECTION  4 


COLLISION  PROCESSES  IN  THE  ACCELERATION  REGION 


4.1  SPECTROMETRIC  EVIDENCE  FOR  COLLISIONS  IN  EXHAUST  BEAM 

4.1.1  FRAO'lENT  IONS 

Evidence  for  collisions  in  the  colloid  thruster  beam  can  be  Inferred 
from  a comparison  of  the  data  available  from  the  electron  impact  spectra 
(El)  for  glycerol  and  the  electrohydrodynamic  spectra  (EH)  using  a Nal- 
glycerol  working  fluid.  El  spectra  are  typically  generated  by  intro- 
ducing vaporized  sample  molecules  into  an  ion  source  chamber  which  are 
bombarded  by  a beam  of  low  energy  electrons.  Ions  produced  from  the 
collisions  are  focused  into  a mass  spectrometer.  Mass  spectra  indi- 
cating the  relative  amounts  and  masses  of  positive  ion  fragments  formed 
by  electron  bombardment  of  glycerol  molecules  have  been  measured. 

EH  spectra  are  measured  by  the  direct  coupling  of  an  electrohydrodynamic 
source  of  charged  particles  (e.g.,  colloid  needle  emitter)  to  a mass 
spectrometer  for  the  detection  of  ionic  species.  This  technique  was 

recently  used  to  study  the  fragment  ions  produced  when  Nal-glycerol  is 

18 

subjected  to  intense  fields  at  a needle  emitter.  Table  III  lists 
some  glycerol  fragment  ions  observed  in  common  by  both  EH  and  El 
spectra  analyses. 

Other  fragment  ions  were  observed,  but  those  listed  in  Table  III  were 
easily  identifiable  and  represent  those  species  having  the  largest  rela- 
tive Intensities  for  both  spectra  techniques. 
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Fragment  ions  with  mass  units  less  than  92  measured  by  EH  spectra  are 
most  probably  the  result  of  the  bombardment  of  neutral  glycerol  mole- 
cules. This  assumption  appears  plausible,  since  the  major  glycerol 
fragment  ions  exhibiting  the  highest  relative  intensities  show  up  in 
both  EH  and  El  spectra.  The  fragment  ions  produced  in  the  El  measure- 
ments are  clearly  formed  by  collisions  of  low  energy  electrons  with 
glycerol  vapor.  The  bombarding  particle  hypothesized  for  the  EH  process 
which  interacts  with  neutral  glycerol  vapor  to  produce  fragment  ions 
could  be  electrons  produced  in  the  accelerating  region.  These  fragment 
ions  could  also  be  produced  by  collisions  of  energetic,  positively 
charged  molecular  ions  with  neutral  glycerol  molecules.  The  probabil- 
ity of  producing  fragment  ions  (m  < 92)  from  collision  between  charged 
droplets  and  neutral  or  molecular  ions  is  assumed  smaller  than  bom- 
bardment of  neutral  molecules. 


In  summary,  evidence  that  fragment  ions  are  produced  in  collisions  and 
not  field  generated  from  emission  sites  is  inferred  from  the  similari- 
ties of  fragments  observed  in  the  EH  and  El  collision  spectra.  Support- 
ing evidence  is  presented  in  the  next  subsection,  where  field  strengths 
necessary  to  generate  various  ionic  species  by  EH  methods  are  examined. 
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4.1.2  FIELD  STRENGTHS  FOR  EVAPORATION  OF  IONS  IN  EHD  PROCESS 


In  the  absence  of  an  electric  field,  the  energy  required  to  evaporate 
an  ion  is  given  by 

+ I - 0 . (22) 

where  is  the  vaporizing  energy,  I is  the  ionization  energy,  and 

0 is  the  work  function  or  the  energy  gained  when  an  electron  is  returned 

to  the  material.  The  relation  between  the  applied  electric  field 

strength  and  the  energy  required  for  field  evaporation  of  singly  charged 

19 

ions  is  expressed  by  the  equation; 

3 1/2 

Q = - (e^  , (23) 

where  the  units  are  e.s.u.  and  e is  the  electronic  charge  and  the  right- 
hand  terms  shows  the  reduction  in  the  total  energy  (Q)  by  the  electric 
field  required  for  the  ions  to  escape.  Tlie  applied  accelerating  field, 
E^,  must  overcome  the  Schottky  mirror  image  force  acting  on  the  molecu- 
lar ion  as  it  escapes  from  the  liquid  substance.  The  ion  escapes  when 
the  energy  barrier,  Q,  is  reduced  to  zero  or  is  very  small  compared  with 
other  terms  in  I'q.  (23)  so  that 

Ev  = -\  (H^  + I - 0)^.  (24) 

e 

Polarization  terms  are  neglected  in  Eq.  (24)  as  they  are  assumed  to  be 
very  small.  If  the  values  of  H^,  I,  and  0 are  expressed  in  electron 
volts,  Eq.  (24),  the  evaporation  fields  in  volts  per  centimeter  can  be 
calculated  from: 

E^  = 6.95  X 10^  (11^  + 1-0)^  (V/cm)  . (25) 

Depending  on  the  value  of  the  applied  field,  E^,  molecular  ions  may  be 
evaporated  in  various  forms. 


In  the  KUD  spraying  process,  ions  of  the  form  C • Na^  (n  is  an  integer), 

+ + '^  + 

G , and  molecular  fragment  ions  (e.g.,  CH2OH  ) tlenoted  by  G^,  have  been 

detected  by  means  of  time-of-flight  or  Ell  spectra  techniques  (sub- 
section 3,2).  If  is  expressed  by  the  appropriate  physical  quanti- 
ties for  ion  formation,  the  evaporation  fields  for  the  ionic  types 
listed  above  can  be  estimated.  The  following  forms  for  are  suggested: 


G • Na"^ 
n 

Ql 

= H 

V 

Type  I 

G+ 

^2 

= H +1-0 

V 

Type  11 

(26) 

Cf 

^3 

- 

Type  e 

where  E is  binding  energy  between  atoms,  and  I,-  is  the  ionization  poten- 
tial  for  a fragment  ion.  For  Type  I ions,  the  term  (1  - 0)  has  been 
dropped,  since  the  accelerating  field  does  not  have  to  provide  the  energy 
for  ionization.  The  Na^  ion  has  been  made  available  previously  via  the 
electrochemical  solvation  process.  Tlius,  the  energy  provided  in  this 
case  to  evaporate  Type  I ions  must  be  sufficient  to  overcome  only  the 
cohesive  or  Van  der  Waal  forces  which  are  represented  by  the  vaporiza- 
tion energy.  The  production  of  a Type  II  ion  under  the  action  of 
intense  fields  is  most  similar  to  Muller's  field  evaporation  of  metal- 
lic ions.  Energy  is  required  not  only  for  the  breaking  of  cohesive 
bonds  arising  from  extramolccular  forces  but  must  also  be  sufficient  to 
ionize  the  glycerol  molecule.  Type  III  ions  or  fragment  ions  differ 
from  Type  II  ions,  since  energy  must  be  supplied  to  fragment  tlie  mole- 
cule by  overcoming  intermolecular  forces.  Intermolecular  bonding  in 
the  glycerol  molecule  is  covalent  in  nature  and  the  bond  dissociation 
energy  for  C-C  is  on  the  order  of  5 to  6 eV.  This  is  10  to  50  times 
strongerthan  dispersive  forces  connecting  molecules  to  each  other. 
Therefore,  the  evaporation  energy,  H^,  for  Types  1 and  I I ions  is 
replaced  by  E„,  the  eond  dissociation  energy,  in  the  activation  energy 

D 

expression.  The  ionization  energy  for  the  glycerol  molecule  is  re- 
placed by  Ij,  the  ionization  energy  for  the  fragment  ion.  I^  will 
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probably  be  somewhat  larger  than  I,  but  not  significantly  and  will  not 
exceed  the  ionization  energy  for  tl>e  oxygen  atom  (13,6  eV). 

The  electric  field  strengths  to  generate  the  various  ions  were  computed 
using  Eq.  (25)  after  substitution  of  the  appropriate  term  previously 
discussed  and  are  listed  in  table  IV  along  with  their  likely  origins. 


TABLE  IV 

FIELD  STRENGTHS  REQUIRED  TO  EVAPORATE  VARIOUS  IONS  BY  EHD  EMISSION 


Type 

H^(eV) 

I,  If(eV) 

Eg(eV) 

E (V/cm) 
V 

Origin 

G • 

n 

Na"^  (I) 

0.9 

— 

-- 

-10^ 

EHD  Generated 

(II) 

0.9 

10 

— 

«.10® 

Collision  Product 

(III) 

-- 

10 

5 

o^ 

0 

1 

Collision  Product 

For  metallic  systems,  0 represents  the  absolute  work  function  that  would 
be  the  energy  gained  by  a metal  in  the  production  of  ions  of  Type  II. 

For  the  condition  of  Type  II  and  IVpe  III  ion  production  in  glycerol 
by  EHD  emission,  energy  should  also  be  gained  by  the  capture  of  an 
electron  into  the  material  following  ionization.  A value  of  10  eV  was 
estimated  for  the  ionization  energy  of  glycerol  for  the  above  calcula- 
tions. An  examination  of  tables  containing  ionization  potentials  for 

organic  molecules  indicates  that  a value  of  10  eV  is  typical  for  organic 

20 

molecules  over  a wide  range  of  structural  complexity. 

The  electric  field  strengths  existing  at  the  tip  of  the  emitters  used 
in  the  TRW  thruster  are  probably  on  the  order  of  10^  V/cm.  In  order 
to  produce  Ions  of  the  type  G^Na^  requires  field  strengths  an  order  of 
magnitude  greater  than  the  fields  at  the  emitter.  Local  field  strengths 
of  these  magnitudes  can  easily  be  achieved  in  the  vicinity  of  emission 
sites  due  to  liquid  spikes.  In  view  of  the  extreme  field  strengths 
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2 3 + 

(10  to  10  higher  over  that  for  the  production  of  • Na  ions) 

required  to  generate  or  glycerol  fragment  ions,  it  is  assumed  that 
these  ions  are  not  generated  by  the  ElID  emission  process.  The  conclu- 
sion is  that  ions  of  Types  II  and  III  observed  in  EH  spectra  are  the 
result  of  collisions  within  the  beam.  It  was  previously  mentioned  that 
glycerol  fragment  ions  were  observed  in  EH  spectra,  but  G^  ions  have 
also  been  observed. 


The  following  system  of  reactions  is  identified  as  collision  processes 
in  the  colloid  thruster  exhaust  beam  that  deserve  attention: 


a.  Collisions  between  molecular  ions  and  charged  droplets 

b.  Collisions  between  molecular  ions  and  neutral  glycerol  molecules 

c.  Collisions  between  molecular  ions  and  electrons 

d.  Collisions  between  charged  droplets  and  electrons 

e.  Collisions  between  neutral  glycerol  molecules  and  electrons 

f.  Charged  droplet  collisions  with  charged  droplets 

Collision  rates  and  energies  involved  in  process  (a)  appear  more  favor- 
able as  a source  of  negative  ion  formation.  The  other  processes  have 
been  examined  and  found  less  likely  as  a source  for  negative  ions  for 
the  following  reasons.  Process  (b),  collisions  between  molecular  ions 
and  glycerol  molecules,  has  a low  probability  of  occurring  because  of 
the  low  molecular  density  only  a few  emitter  radii  away  from  the  emitter 
surface.  In  addition,  neither  of  the  species  contains  the  iodine  needed 
to  produce  the  negative  ion.  Process  (c),  collisions  between  molecular 
ions  and  electrons,  suffers  from  the  same  conditions  of  process  (b). 
Process  (d),  collision  between  charged  droplets  and  electrons,  will  be 
discussed  in  the  next  section.  Process  (e),  collisions  between  mole- 
cules and  electrons,  also  lacks  the  source  of  negative  ions.  Process  (f), 
collisions  between  charged  droplets,  could  possibly  produce  negative 
ions,  depending  upon  the  collision  energy.  This  process  involves  a 
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a range  of  relative  energies  related  to  the  range  of  droplet  charge-to- 
mass  ratios,  but  this  is  encompassed  by  process  (a),  where  molecular 
ions  have  the  highest  charge -to -mass  ratios  in  the  distribution  and 
the  rest  of  the  distribution  is  represented  by  a single  particle. 

When  the  entire  distribution  is  taken  into  account,  the  results  are 
found  not  to  differ  greatly  from  the  mean  droplet  model.  Collision 
between  droplets  can  occur  and  can  produce  negative  ions  according  to 
the  analyses  discussed  in  the  following  sections.  This  will  simply 
add  to  the  results  found  and  does  not  invalidate  the  model. 

The  background  development  of  droplet  energetics  discussed  in  sub- 
section 3,4  aids  in  understanding  charged  droplet  breakup  or  evapora- 
tion when  energy  is  dumped  into  a droplet  via  a collision  with  a molec- 
ular ion.  This  section  will  treat  aspects  of  collisions  involving 
droplets,  since  droplet  breakup  suggests  several  possible  mechanisms 
for  the  production  of  negative  ions.  An  analysis  of  electron  colli- 
sions with  droplets  was  undertaken  on  the  assumption  that  some  electrons 
reach  the  accelerating  region  either  from  the  chamber,  from  thruster 
electrodes,  or  that  they  are  formed  within  the  region  through  beam 
collisions . 

A discussion  of  electron  collisions  is  given  in  the  next  subsection, 
followed  by  an  analysis  of  molecular  ion  collisions  with  charged  droplets. 


4.2  ELECTRON  COLLISIONS  WITH  CHARGED  DROPLETS 


An  empirical  relationship  relating  the  range  of  monoenergetic  electrons 
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in  materials  as  a function  of  energy  is  given  by  the  expression: 

R = aE^"  (mg/cm^)  , (27) 

where  a and  n are  empirical  constants  and  E^  is  the  energy  of  the  elec- 
trons in  MeV.  The  constants  have  the  following  values  for  the  energy 
indicated ; 
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a 


710,  n 


= 1.72  (0.001  ^ 0.3  MeV) 


The  range  in  Eq.  (27)  is  in  units  of  mass  per  unit  area,  and  conversion 
to  penetration  distance  in  units  of  centimeters  is  accomplished  through 
the  relationship  R = px  to  yield: 


a X,  V 

X = — E (cm) 

P o 


where  pis  the  density  of  the  absorbing  material.  For  glycerol,  the 
penetration  distance  for  electrons  as  a function  of  energy  becomes: 

X = 0.563  E . (29) 

o 

Figure  10  is  a plot  of  the  penetration  (x)  in  glycerol  as  a function  of 
electron  energy  in  kV.  < 


For  comparison  purposes,  the  electron  range  equivalent  to  the  diameter 
o _6 

of  50A  radius  (D  = 10  cm)  droplet  is  sho^im  on  the  plot  and  occurs 

for  an  electron  energy  of  450  to  500  eV.  Therefore,  electrons  col- 

o . ^ 

liding  with  a 50A  droplet  with  energies  greater  than  0.5  kV,  will  pass 

through  the  droplet.  For  example,  a 1 kV  electron  would  penetrate 

glycerol  a distance  of  4 x 10  ^ cm,  which  is  four  droplet  diameters. 

For  energies  less  than  0.5  kV,  electrons  may  be  captured  by  the  droplet, 

depending  upon  the  angle  of  incidence  of  the  collision.  For  electrons 

that  originate  in  the  chamber  downstream  from  the  thruster  and  reach  the 

accelerating  region,  capture  would  take  place  only  in  the  outer  fringes 

of  the  accelerating  region. 


The  energy  loss  per  unit  path  length  due  to  ionization  for  electrons 
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moving  through  a glycerol  droplet  can  be  estimated  by: 


(£) 


NZ  1 
eft 


crgs/cm, 


t 


f 

1 

I 


where 


= electronic  charge 
= mass  of  electron 
= electron-'Virlocity 

3 

= number  of  atoms/cm 
= effective  value  of  atomic  number 

= geometric  mean  ionization  and  excitation  potential 


Ihe  rate  of  energy  loss  was  calculated  for  the  special  case  wliere  the 
electron  traverses  the  diameter  of  a mean  droplet.  The  results  are 
shown  in  Figure  11,  using  the  following  values  for  glycerol: 
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1.15  X 10  atoms/cm 


= 70.7  eV 


In  Eq.  (30),  is  defined  as  the  geometric  mean  of  the  ionization  and 

excitation  potentials  of  the  material.  The  I value  for  glycerol  has 
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been  determined  experimentally  as  70.7  eV.  Tlie  effective  value  for 
the  atomic  number  of  glycerol  was  calculated  on  the  basis  that  glycerol 

3 

is  a medium  which  is  a mixture  of  N2,  and  atoms/cm  of  carbon, 
hydrogen,  and  oxygen  having  atomic  numbers  Z^,  Z2,  and  Z^,  respectively. 
As  seen  in  Figure  11,  the  maximum  energj'  deposited  in  a droplet  of  5oX 
radius  is  approximately  470  eV,  corresponding  to  the  electron  energy 
just  needed  to  escape  the  droplet.  This  energy  represents  one-sixth 
of  the  total  droplet  binding  energy  and  could  be  quite  effective  in 
fragmenting  the  droplet.  Note  from  Figure  11  that  the  energy  deposited 
in  the  droplet  by  electrons  increases  linearly  with  the  electron  energy 
in  the  region  for  total  caputre  of  the  electron  by  the  droplet.  After 
the  electron  has  enough  energy  to  pass  completely  through  a droplet, 
the  amount  of  energy  deposited  decreases  with  increasing  electron  energy. 
High  energy  electrons  will  not  deposit  sufficient  energy  to  produce 
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many  Ionization  events  within  the  droplet.  Except  for  Ionization  at 
the  surface,  secondary  electrons  produced  In  the  Interior  of  the  drop- 
let probably  lack  sufficient  energy  to  escape  from  the  charged  droplet. 
The  secondary  electrons  produced  by  Ionization  In  a droplet  that  do 
not  escape  will  finally  dissipate  their  energy  In  the  form  of  heat. 
Collisions  between  electrons  and  droplets  can  result  In  the  production 
of  negative  Ions  similar  to  the  model  discussed  In  the  next  section. 

f 

But  the  number  of  electrons  In  the  acceleration  region  Is  very  small 
compared  with  the  other  species  because  any  significant  electron  current 
will  result  In  excessive  emitter  tip  heating  due  to  electron  bombardment 
Electron  trap  electrodes  have  been  used  In  order  to  avoid  this  latter 
problem.  Therefore,  no  significant  negative  Ion  production  Is  expected 
by  electron  bombardment  of  charged  droplets. 

4.3  MOLECULAR  ION  COLLISIONS  WITH  CHARGED  DROPLETS 


The  development  of  a collision  model  that  satisfactorily  predicts  the 

8 9 

production  of  negative  ions  of  the  order  10  to  10  /sec  will  require  an 
examination  of  the  meaningful  collision  processes  in  the  beam.  The 
charged  droplet  is  a carrier  of  iodine  ions  and  iodine  atoms,  which 
have  a large  electron  affinity  (3.06  eV).  This,  coupled  with  the 
instability  of  the  droplet  for  collislonal  energies  exceeding  several 
hundred  eV,  leads  to  the  selection  of  the  interaction  between  charged 
droplets  and  molecular  ions  as  the  most  critical  collision  process  for 
study.  It  is  convenient  to  examine  the  collision  between  ion  and  charged 
droplets  In  terms  of  the  collision  rate  given  by: 


I = I p ax  (collision/sec)  , (31) 

c o 

where  I Is  the  molecular  Ion  current  (ions/sec),  pis  the  density  of 

° 3 2 

charged  droplets  (#/cm  ),  o is  the  collision  cross  section  (cm  ),  and 

X is  the  interaction  distance  (cm).  The  cross  section  for  this  collision 
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is  taken  as  the  geometric  cross  section  of  the  droplet  which,  in  general, 
will  be  a function  of  tlie  number  of  molecules  in  a droplet  and  is  char- 
acterized by  the  symbol  o^.  The  justification  for  using  the  geometric 
cross  section  was  based  on  an  analysis  of  charged  particle  scattering. 

Consider  the  case  for  a positively  charged  ion  approaching  a positively 

charged  droplet.  It  will  experience  a coulomb  repulsive  force  referred 

to  as  Rutherford  scattering.  This  could  deflect  the  ion  sufficiently 

to  prevent  a direct  collision  under  certain  conditions.  In  effect,  it 

would  decrease  the  collision  cross  section  below  the  geometric  cross 
2 

section  (nR  , where  R is  the  droplet  radius). 

Calculations  were  made  using  Rutherford  scattering  theory  to  determine 
the  energy  at  which  ions  are  sufficiently  deflected  to  avoid  a direct 
collision  with  the  mean  specific  charge  droplet.  At  high  enough  ener- 
gies, ions  collide  with  the  droplet  as  long  as  the  impact  parameter  is 
equal  to  R or  less.  Analysis  shows  that  the  assumption  for  the  geo- 
metric cross  section  is  valid  for  ion  energies  down  to  approximately 
100  eV.  The  collision  cross  section  drops  to  one-half  the  geometric 
value  near  20  eV.  Since  most  of  the  critiOal  collisions  will  occur  at 
energies  well  above  100  eV,  it  is  sufficient  to  use  the  geometric  cross 
section  for  collisions  between  ions  and  droplets. 

As  the  charged  droplet  density  is  not  constant  in  the  accelerating 
region,  the  droplet  density  function  0(x)  derived  in  subsection  3.3 
will  be  used  in  the  expression  for  collision  rates.  Thus, 

Ic  = / x^  (32) 

which  becomes,  after  substitution  of  the  density  function  Eq.  (12), 


1 0„ 
o g 


‘1  (2c  [v^  - v(x)]3^ 


This  expression,  therefore,  describes  the  total  number  of  collisions 
that  occur  over  the  interaction  distance  X2  ” Since  the  integrand 

of  Eq.  (33)  contains  V(x),  the  limits  of  integration  can  be  converted 
from  distance  to  potential.  In  terms  of  collision  energet^s,  it  is 
useful  to  know  the  total  number  of  collisions  that  occur  pe^r  unit  time 
after  the  ion  has  been  accelerated  through  a given  potential  difference. 
This  was  accomplished  bv  solving  the  potential  function  foirx  in  Eq.  (1) 


This  was  accomplished  by  solving  the  potential  function  fo^: 
of  Section  2,  With  the  appropriate  substit^ition^o^flfe^e^jt 
the  collision  rate  can  be  reduced  to:  ^ 


ssion  for 


- 


dy, 


where  K is  a constant  and  the  integral,  known  as  Dawson's  integral,  is 
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tabulated  in  handbooks.  Tlie  constant  K is  given  by: 

'“©r ' 


K “ r I 

o o g Z L V 


Collision  rates  were  calculated  for  the  following  set  of  conditions: 


V = 13  kV 

o 


C = 10  C/kg 


R = 1.57  X 10 


r = 7.8  X 10 

o 

12 

I = 6 X 10  ions/sec 

o 

19  2 

0 = 2.04  X 10  droplcts/m  -sec 

= 7.85  X lO"^ 
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The  results  are  shown  in  Figure  12,  where  the  total  number  of  molecular 
lon/droplet  collisions,  wtiich  occur  after  the  ion  lias  been  accel- 
erated through  a given  potential  is  plotted  as  a function  of  that 
potential.  For  example,  molecular  ions  that  have  been  accelerated  from 

g 

the  emitter  through  a potential  of  3 kV  will  then  undergo  7.6  x 10 
collisions  with  droplets  before  leaving  the  accelerating  region.  Note 
that  for  high  energy  collisions  involving  molecular  ions  that  have  been 
accelerated  through  at  least  10  kV,  the  collision  rates  are  still  on 

g 

the  order  of  10  /second.  The  collislonal  energy  for  the  molecular  ion/ 
droplet  system  is  approximated  by 


E = 1/2  U vt  , 

c ion 


V.  » , 

ion  D 


(36) 


where  is  the  reduced  mass  and  v.  , the  velocity  of  the  molecular  ion, 

ion 

is  much  greater  than  the  droplet  velocity  v^^.  Further,  since  the  drop- 
let mass  is  much  greater  than  the  mass  of  the  molecular  ion,  M.  , 

ion 

Eq.  (36)  becomes; 


E = 1/2  M,  V,  = qV. 

c ion  ion 


Thus,  the  collisioral  energy  in  electron  volts  involving  singly  charged 

ions  and  heavy  charged  droplets  can  be  estimated  by  the  potential 

through  which  the  ion  has  been  accelerated  at  the  time  of  the  collision. 

Collisions  between  droplets  and  ions  that  have  been  accelerated  through 

10  kV  involve  collision  energies  three  times  greather  than  the  droplet 

total  binding  energy  (E  ~ 3E  ),  From  Figure  12,  the  number  of  colli- 

^ ^ 8 

sions  at  energies  of  E^.  or  greater  is  larger  than  10  . Assumed  produc- 
tion rates  of  1 negative  ion  per  collision  would  then  erode  about 
1 mil/year  at  sputtering  yields  even  as  low  as  1 atom/ion. 
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4.4  HVAPOKATU'N  KKOM  CllARCKI)  OKOIM.KTS 

The  char{;od  liquiil  droplots  fU’iior.itoJ  nml  nccploratod  tliroiif’li  tlit> 
Intense  electric  fieltl  at  the  emitter  can  uncicrjjo  various  interactions 
with  the  fields  and  other  particles.  Kacli  of  tlieso  interactions 
affects  the  stability  of  tlie  droplets  and  may  lead  to  particle  breakup. 
Particle  breakup  could  also  occur  due  to  the  inherent  instabilities 
initiated  by  conditions  other  than  collisions  and  interactions. 

Wlien  cliarged  droplets  arc  generated,  they  are  assumed  limited  by  tlie 
Rayleigh  criterion.  Tliat  is,  the  droplet  siac  and  total  ciiarge  are 
limited  so  that  the  maximum  charge-to-mass  ratio  for  glycerine  is 
given  by: 


C = 3.54  X lO"^  , (37) 

wl»ere  C is  i n units  of  coulombs  per  kilogram  and  r is  t!ie  droplet  radius 
in  meters.  Prop  lets  can  he  generated  witli  cliarge-to-mass  ratios  equal 
to  or  less  than  the  value  of  C in  liq . (37).  However,  since  the  droplets 
undetgo  oscillations  upon  generation,  they  will  sliov^  some  instability 
because  the  Rayleigh  criterion  is  based  upon  tlie  droplets  being  splu'ri- 
cal.  Tims,  tlie  temleney  toward  droplet  breakup  will,  in  part,  depend 
upon  the  oscillatory  behavior.  The  oscillation  frequency  is  related 
to  viscosity  which,  in  turn,  depends  strong  upon  temperature.  Unfor- 
tunately, little  is  presently  known  about  microscopic  droplet  oscilla- 
tions and  the  cffect:s  of  oscillations  upon  stability  during  the  genera- 
tion and  acceleration  processes.  Similarly,  little  is  known  about  the 
droplet  temperature  upon  generation.  It  is  suspected  that  the  tempera- 
ture is  far  greater  than  the  emitter  tip  temperature  (typically  at  25''c) 
as  a result  of  the  current  flow  during  the  generation  process. 

Neutral  molecule  evaporation  from  the  chargeil  droplet  will  occur  at  a 
rate  dependent  upon  the  temperature.  This  evaporation  could  cause  a 
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rapid  decrease  in  the  droplet  size  so  that  the  radius  decreases  and 
the  charge-to-mass  ratio  increases,  resulting  in  a departure  from  the 


liniiatioiis  of  the  Rayleigh  criterion  as  expressed  by  Eq.  (37).  This 

25 

phenomenon  was  described  by  Abbas  and  Latham  in  which  evaporation 


and  successive  droplet  instabilities  were  observed. 


As  discussed  previously,  the  droplet  temperature  can  undergo  a drastic 
increase  as  a result  of  a collision  with  an  energetic  ion.  This  could 
cause -partial  or  even  total  evaporation.  If  the  evaporation  time  is 
comparable  to  the  time  spent  in  the  accelration  region,  the  probability 
of  generating  negative  ions  which  can  return  to  the  emitter  is  greatly 


enhanced.  The  time  required  for  heat  to  transfer  throughout  the  droplet 

r‘  f .11  1 . ./.**11  1 _ l _ 


from  a collision  would  be  about  10  seconds,  assuming  that  heat  trans- 
fer propagates  at  the  velocity  of  sound  for  such  a small  droplet.  This 


is  a short  time  comparted  to  the  acceleration  region  transit  time 


(~  10  sec).  It  would  be  of  interest  to  compare  tlie  evaporation  time 
for  various  temperatures  with  the  above  transit  time. 


The  time  for  the  droplet  to  evaporate  90  percent  of  ots  mass  was  cal- 
culated and  defined  as  the  evaporation  time  (i).  This  time  decreases 
with  increasing  temperature  in  accordance  with  the  evaporation  rate. 

The  calculation  was  simplified,  since  it  is  for  comparative  purposes, 
and  some  factors  were  neglected.  The  temperature  is  lowered  by  evapora- 
tion which  increases  the  evaporation  time.  Tlie  evaporation  rate  goes 
up  as  the  droplet  size  decreases,  which  decreases  the  evaporation  time. 
These  two  important  parameters  will  be  discussed  later  and  their  effect 
upon  the  calculations  considered. 


The /change  per  unit  time  in  the  number  of  molecules  (n)  in  the  droplet 


due  to  evaporation  is  given  by 


tt  = ■V(^>  ' 
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where  is  the  molecular  evaporation  rate  in  molecules/m  -sec,  and 


* 

A(t)  is  the  droplet  area  which  varies  with  time.  Assume  the  droplet 


radius  is  R = 50A,  and  the  volume  of  a molecule  in  the  droplet  is 
o 3 

V = (5A)  , Kq.  (38)  becomes 


dn 


2/3  " dt  . 


(39) 


This  equation  is  integrated  from  n , the  initial  number  of  molecules, 


to  0.1  n on  the  left,  and  from  0 to  t on  the  right.  After  integration, 
o 


the  evaporation  time  is  found  to  be 

1/3 


3(0.A64)  n 


(40) 


Considering  droplets  raised  to  various  temperatures  at  generation,  due 


to  collisions  or 

for  any  other  condition,  the 

evaporation  times  and  the 

evaporation  rates 

are  given  in  Table  V. 

1 

TABLE  V 

• 

EVAPORATION  TIME 

Tempera  l ure 
^cc Icius 

Evaporation  Rale 
J^(mol/m^-sec) 

Evaporation  Time 

T (sec) 

20 

lo^o 

• 2.09  X lO"’ 

50 

21 

1.8  X 10 

CM 

1 

O 

X 

i 

100 

23 

1.4  X 10 

-4  • 

1.49  X 10 

. 1 

150 

24 

3.09  X 10 

6.8  X lO"^ 

! 

200 

25 

2.7  X 10 

7.7  X 10“^ 

1 

1 

' 290 

4.46  X 10^^ 

-ft 

4.68  X 10 

From  the  results  given  in  Table  V,  illustrated  in  Figure  13,  droplets 
with  temperatures  well  above  200°C  will  undergo  total  evaporation  within 
the  acceleration  region.  But  partial  evaporation  is  possible,  resulting 
in  Rayleigh  instabilities  down  to  temperatures  close  to  100°C,  since 
approximately  30  percent  of  the  droplet  will  evaporate  in  one-fifth  of 
the  calculated  evaporation  time.  In  theory,  collisions  with  ions  having 
energies  of  a couple  of  kV  or  more  will  undoubtedly  result  in  droplet 
disintegration  within  the  acceleration  region.  Below  these  energies, 
droplet  fragmentation  will  readily  occur  with  perhaps  loss  probability 
of  negative  ion  production. 

Since  the  decrease  in  temperature  due  to  evaporation  was  not  included, 

the  evaporation  times  would  be  increased  above  the  calculated  values. 

To  Include  cooling  effects,  the  heat  of  vaporization,  the  specific  heat, 

and  the  thermal  evaporation  rate  must  be  carried  as  > '.riables,  and 

these  functions  do  not  have  simple  analytic  expressions.  However,  the 

increase  in  evaporation  time  is  somewhat  balanced  by  a corresponding 

2 6 

increase  in  evaporation  rate  with  decreasing  droplet  size.  This  is 
due  to  the  droplet  curvature  approaching  molecular  dimensions  and  loss 
of  bonds  holding  the  molecules  on  the  surface.  Iv/lien  tlie  droplet  radius 
decreases  to  one-fifth  of  its  initial  size,  the  evaporation  rate  in- 
creases by  close  to  a factor  of  three. 
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EVAPORATION  TIME  (sec) 


■n.-JB®' 
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SECTION  5 

EFFECT  OF  A DISTRI BimON  OF  DROPLETS  ON  THE  COLLISION  RATE 


The  charged  glycerine  particles  emitted  from  the  EHD  source  can  be 

classified  into  two  separate  groups.  The  first  group,  which  consti- 

tures  approximately  91  percent  of  the  total  number  of  charged  particles, 

consists  of  charged  liquid  droplets.  The  second  group  consists  of 

singly  charged  molecular  ions.  In  the  first  group,  the  mean  value  of 

4 

the  charge-to-mass  ratio  (q/m)  is,  q/m  = 10  C/kg,  and  in  the  second 
group  it  is,  q/m  = 4.62  x 10^  C/kg.  Since  both  groups  are  accelerated 
by  the  same  electrical  potential,  the  velocity  of  the  tons  will  be 
approximately  seven  times  greater  than  that  of  the  dr-^plets.  The 
faster  moving  ions  constitute  an  electric  current  which  passes  through 
the  aggregate  of  slower  moving  droplets.  It  is  of  interest  to  deter- 
mine the  collision  rate  of  the  ions  as  a function  of  the  distance 
traversed  by  them  through  the  relatively  motionless  droplets. 


5 . 1 DEFINITION  OF  ION  COLLISION  RATE 


Let  I^  denote  the  initial  magnitude  of  the  molecular  ion  current,  and 

let  I denote  the  current  after  it  has  moved  a distance,  x,  through  the 

"fixed"  group  of  droplets.  Now,  the  number  of  singly  charged  molecular 

ions  leaving  the  ion  source  in  time,  At,  is,  (I  /e)At,  where,  e = charge 

19  ° 

on  electron  = 1.6  x 10  coulombs.  The  niunber  of  singly  charged  molec- 
ular ions  that  survive  a distance  x (do  not  make  collisions)  through 
the  droplets  in  time  At  is,  (I/e)At.  Ihe  number  of  ions  that  collide 
with  the  droplets  is  then,  (I  /e)At  - (I/e) At  and  the  collision  rate, 

P,  is  this  quantity  divided  by  At.  Thus, 

I - l(x) 

P(x)  = , collisions/sec.  (41) 
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To  determine  0(x) , an  expression  Cor  I(x)  must  first  be  derived.  This 
is  done  in  the  following  paragraphs. 


5.2  MOLECULAR  ION  CURRENT 


5.2.1  DROPLET  VELOCITY 


It  is  assumed  that  the  liquid  droplets  are  spherical  and  that  their 
radii  are  generally  unequal.  Let  a(r^)  denote  the  collision  cross  sec- 
tion of  those  droplets  that  have  a radius  r^,  and  let  n(r^)  denote  the 
particle  number  density  of  these  same  droplets.  Then,  the  reduction  in 
the  Ion  current,  dl,  after  traversing  a distance,  dx,  is  given  by 


dl  = -I  j^^(r^)n(r^)j  dx  , 


where  the  quantity,  !E]CT(r^)n(r^),  represents  the  "absorption  coefficient" 
of  the  tons.  Here, 


r(r)  = ■nr'^ 


n(r)  = 


where  A = total  cross-sectional  area  of  the  ion  beam,  v = velocity  of 

. ^ 

those  droplets  having  a radius  r,  and  = flow  rate  (particles/sec) 
of  those  droplets  having  a radius  r. 

The  molecular  ions  and  charged  liquid  droplets  are  both  accelerated 
through  a potential  difference,  V^-V(x),  where  = potential  at  the 
charged  particle  emitter,  and  V(x)  = potential  at  a distance  x from 
the  emitter  along  the  axis  of  symmetry,  A good  approximation  of  this 
potential  taken  from  Section  3 is  given  by 
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V(x) 


(X  < a)  , 


(45) 


o -in  (a/b)  ’ 


-3  -4 

where,  a = 1.57  x 10  meters,  b = 3.9  x 10  meters,  and,  = 13,000 
volts.  Note  that,  V(0)  = V^.  ITie  velocity  acquired  by  those  droplets 
having  a charge-to-mass  ratio,  q/m,  after  being  accelerated  through  a 
distance  x is 

_ r T 1/2 

V = J2  1_V^  - V(x)J  (q/m)  . (46) 


It  is  experimentally  verifiable  that  q/m  is  related  to  the  droplet 
radius,  r,  fairly  accurately,  by 
.1/2 

. (47) 


a = 

m 


3(e  yY 
-O. 


3/2 


or 


where,  Cq  ® permittivity  of  free  space  = 8.854  x 10  farads/meter, 
V = surface  tension  of  droplet  = 1.3  ne\rtons/meter , P = density  of 

3 

droplet  = 8000  kg/meter  . Combining  Eqs.  (46)  and  (47)  yields  the 
desired  droplet  velocity. 


(48) 


Since  v^  = 0 when  x = 0,  therefore,  from  Eq.  (44)  the  particle  number 
density  becomes  infinite  at  x = 0,  Since  this  is  not  a realistic 
result,  the  minimum  allowable  value  of  x,  denoted  by  x^,  must  be  made 
greater  than  zero.  A practical  value  of  x^  is  assruned  to  be  that  which 
satisfies  the  equation. 


V - V(x  ) = 100  volts  , 

o o 


(49) 


54 


since  the  faster  moving  ions  only  begin  to  break  up  the  droplets  into 
separate  ions  (by  collision)  after  the  molecular  ions  have  been  accel- 
erated through  a potential  difference  of  approximately  100  volts.  Sub- 
stituting Eq.  (45)  into  Eq.  (49)  yields 

X = 4.23  X 10  ^ meters, 

o 

5.2.2  DROPLET  FLOW  RATE 


Experimental  results  suggest  that  the  charge-to-mass  ratio,  q/m,  of  the 
droplets  is  normally  distributed.  Thus,  if  f(q/m)  denotes  the  proba- 
bility density  function  of  q/m,  then 


f (q/m) 


(50) 


where,  S and  q/m  are  the  standard  deviation  and  mean,  respectively,  of 
the  distribution.  As  stated  earlier,  q/ra  = 10^  C/kg.  Also,  S ^ 3000  C/kg. 
Since  the  droplet  radius,  r,  is  related  to  q/m  through  Eq.  (47), there- 
fore, the  probability  density  function  of  r,  denoted  by  g(r),  can  also 
be  obtained.  It  is  given  by 


g(r)  = f(q/m) 


d (q/m) 


dr 

1/2 


-5/2 


S(2q) 


1/2 


r 

exp  < - — j- 
2S 


'9e  Y 
o 


1/2 


-3/2 


(51) 
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Note  that  g(0)  =0.  If  N denotes  the  total  flow  rate  of  the  droplets, 
then  it  follows  that. 


r+Ar 


Nj.  = N J g(u)  du  = Ng(r)Ar 


(52) 


where,  Ar  is  an  infinitesimal. 

Since  N represents  the  total  flow  rate  of  the  droplets,  therefore,  the 
electric  current,  i,  generated  by  the  motion  of  the  charged  droplets 
can  be  expressed  by 

1 = N q 


where  q denotes  the  mean  charge  of  the  droplets.  Now  the  total  current, 
I^,  leaving  the  emitter  is  given  by 


where  is  the  initial  molecular  ion  current.  Since  i is  approxi- 
mately ten  times  the  magnitude  of  I^,  therefore,  i = 10  1^/11.  Let, 
I^/A  = J^,  denote  the  total  current  density  at  the  ion  source.  Hence, 


N ^ 10  ^ ;Jt 

A 11  ^ ^ 


(53) 


The  quantity  can  be  measured.  The  mean  droplet  charge,  q,  is  deter- 
mined as  follows.  Since  the  droplets  are  assumed  to  be  spherical,  the 
droplet  mass,  m,  is  then  given  by,  m = 4rrr  p/3.  By  combining  this 
equation  with  Eq.  (A7),  it  is  easy  to  show  that 


P q/m 


(54) 
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The  probability  density  function  of  q,  denoted  by  h(q),  is  then  given  by 


h(q) 


f (q/m) 


d(q/m)  ^ 
dq 


12TT^e 

SI 

2 


pq 


f(q/ra)  , 


(55) 


where,  f(q/m)  is  expressed  in  terms  of  q,  via  Eq.  (54).  By  definition, 
the  mean  value  of  q in  the  interval,  qj^  ^q  ^ given  by 


q =•  J qh  (q)dq 


where,  h(q)  is  given  by  Eq.  (55),  and 


(56) 


12TTy'e  , 

o i 

P (q/m)i  ’ 


12TTYe^ 

^2  " p (q/m)2 


The  values  assigned  to  (q/m)^^,  and  (q/m)^  should  be  chosen  so  that  they 
include  practically  all  possible  values  of  q/ra.  Since  q/m  is  normally 
distributed,  therefore,  the  following  assigned  values  of  (q/m)j^  and 
(q/m)2  will  bracket  99.73  percent  of  all  possible  values  of  q/m: 


(q/m)j^  = qTm  + 3S  = 19,000  coulombs/kg 

(57) 

(q/ra)2  = q/m  - 3S  = 1,000  coulombs/kg 


By  introducing  a change  of  variable  from  Eq.  (54),  and  making  use  of 
Eq.  (55),  it  is  easy  to  show  that  Eq.  (56)  becomes. 


19,000  i2TtY€  , 

q =1'  f(q/m)d(q/m)  , (58) 

1,000 

where,  f(q/m)  is  given  by  Eq.  (50). 
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5.2.3  MOLECULAR  ION  CURRENT 


By  combining  Eqs.  (53),  (52),  (A8) , (44),  and  (43),  one  obtains  the 
following; 

2\^'  - 


^ r..  11 

r 4 

V - V(x)  J r g(r)dr  , 
o 


\9Ye 


where,  T2  > and 


■.  ■&)' 


(q/m). 


(q/m) 


2/3  ' 


By  introducing  a change  of  variable  from  Eq.  (47)  and  making  use  of  the 
first  part  of  Eq.  (51),  it  follows  that: 


r2  11 

J*  r^  g(r)dr  = I — ) , J (q/m) 

\ p / 1,000 


11/12  19,000 


f(q/m)d(q/m)  (60) 


Now,  integrating  Eq.  (42)  yields 


I exp 
o 


- J &(r^)n(r^)dx 


Substituting  Eqs.  (60)  and  (59)  into  Eq.  (61)  finally  yields  the  desired 
expression  for  the  molecular  ion  current: 


I(x)  = exp 


2/3  19,000 


fTTjT  /9Ye 

[ — 2“  ) J (9/"') 

? \ p J 1,000 


-11/6 


f(q/m)d(q/m)  . 


(62) 


p dx 

. [V„  - V(x)]^ 


X o 
o 


where,  x^  = 4,23  x 10  meters  and  V(x)  is  given  by  Eq.  (45). 


5.3  ION  COLLISION  RATE 


Recall  that  the  ion  collision  rate  is  given  by  Eq.  <'41),  where  I(x)  is 


given  by  Eq.  (62),  A realistic  value  of  the  total  ion  current  is, 


,"6 


I « 9.6  X 10  amperes,  and  the  cross-sectional  area  of  the  Ion  beam 

^ -3  2 

is,  A = tt(7  X 10  inch)  , Hence, 


J = ~ = 96.67  amperes/meter  , 

1 A 


From  the  values  of  Y>  and  p given  in  subsection  5.2,  it  follows  that 


’9Ye 


2/3 


X 2 A /3  2 

= 1.3785  X 10  (coulombs/kg)  meters  , 


Next,  employing  the  trapezoidal  rule  (numerical  integration),  the 
integrations  needed  to  determine  q,  given  by  Eq.  (58),  and  the  first 
Integral  in  Eq.  (62),  yield  the  following: 
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Employing  the  values  of  a,  b,  given  in  subsection  5.2,  we  have: 

-6  h 

2 bk  =•  8.0732  x 10  (coulombs/kg)^  sec  , 

z(x)  = [1.39268  + 2n(636.94  x +0.2484)]^  , 

z = 0.103697 
o 

Note  that,  since,  x ^ x ^ a,  therefore, 

o 

z s:  z(x)  1,27  . 
o 

Hence,  the  tabulated  values  of  F(y),  for  0 ^ y S ''.00,  given  in  the 
above  reference,  will  be  sufficient  for  evaluating  the  integral  in 
Eq.  (65).  Substituting  all  of  the  foregoing  in  Eq.  ^62)  yields  the 
desired  molecular  ion  current: 


' 2 

■ 

I(x)  = exp 

- 2.61  (lO"^) 

e^  fCz(x)]  - 1.01  F(z  ) 

0 

The  ion  collision  rate  is  then  obtained  by  substituting  Eq.  (66)  into 
Eq.  (41). 

Before  the  collision  rate,  P,  can  be  calculated,  the  initial  ion  current, 

I , must  be  evaluated.  Recall  that  i = 10  I and  i = 10  I_/ll,  so  that 
o o T 

I = I_/ll  = 9.6  X 10  /ll  = 8.727  x 10  amperes.  The  lower  solid 
o T 

curve  in  Figure  14  shows  a plot  of  the  ion  collision  rate  versus  the 

axial  distance,  x,  from  the  tip  of  the  emitter.  Here,  the  collisions 

between  the  ions  and  the  droplets  occur  all  along  the  axial  distance 

from  X to  X,  Note  that  the  total  collision  rate  therefore  occurs  at 
o 

X = a,  i.e.,  at  the  end  of  the  acceleration  region  and  it  is  equal  to 

g 

3.35  X 10  collislons/sec . 
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. DISTANCE  FROM  TIP  OF  ANODE,  X (meters) 

i 


Figure  14.  Molecular  Ion  Collision  Rate  versus  Axial  Distance  from 
tlie  Eraitte’-  Tip 


I 


Recall  that  tfhe  standard '^dpviation  of  the  probability  density  function 
of  q/m,  give\i  by  Eq.  ^^0)  J is  S = 3,000  coulorobs/kg.  What  effect  will 
an  increase  In  S hav^  on  the  ion  collision  rate?  To  answer  this  ques- 
tion, let  S = 3,300  coulombs/kg.  Then  the  upper  and  lower  limits  of 

i ' I 

integration/ in  Eqi.  (58)  ^nd  (62)  become  19,900  coulombs/kg  and 

. / 

100  coulomb/s/kgy  respecti/vely . Hence, 

' / . r 

— I -18 

q = / 6.<4868  x 10!  - coulombs. 


19,960  5 

r I (q/m)  f (d7m)d(q/m)  = 1,2854  x 10  ^ (coulombs/kg) 

, nn  ' ' 


-11/6 


so  that  |Eq.  (62)  then  becomes 


z^(x). 


l(jc)  = exp  -4.;'73(1(?  ^)  ie^  [z(x)]  - 1.01  F(z^) 


Substituting  the  above 'in  Eq.  (41)  yields  the  ion  collision  rate.  The 
middle  solid  curve  in  Figure  14  shows  the  ion  collision  rate  versus  the 
distance,  x.  Note  that  when  x = a,  the  collision  rate  is  81  percent 
greater  than  that  for  the  case  S\  = 3,000  coulombs/kg.  Thus,  a 10  per- 
cent j increase  in  S causes  an  81  percent  Increase  in  the  collision  rate 
at  x!  = a.  \ 


5.4/  APPROXIMATE  ION  COLLISION  RATE 


Al/;hough  Eq.  (66)  gives  an  accurate  description  of  the  molecular  ion 
; \ 
ci^rent,  it  is  worthwhile  to  consider  an  approximate  representation  of 

this  current  as  discussed  in  Section  4.  The  approximation  can  be 

readily  derived  by  making  the  following  two  assumptions: 
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a.  All  the  droplets  have  the  same  charge-to-mass  ratio,  given  by, 
q/m  = 10^  coulombs/kg. 

b.  The  velocity  of  the  droplets  is  independent  of  x,  i.e.,  the 
velocity  is  constant. 

As  a consequence  of  (a),  and  by  virtue  of  Eq.  (47),  all  droplets  will 

-9 

have  the  same  radius,  given  by  r = 2.5295  x 10  meters.  As  a conse- 
quence  of  (b),  Eq.  (46)  becomes,  v = (24V)  (q/m)  , where,  AV  = 13,000  volts, 
so  that,  V = 16,125  meters/sec.  As  a further  consequence  of  assumptions 
(a)  and  (b),  both  the  collision  cross  section  and  the  particle  number 
density  of  the  droplets,  given  by  Eqs.  (43)  and  (44),  respectively,  are 
constant.  Hence,  Eq.  (42)  becomes 

dl  = -landx  , 

which  readily  integrates  to 

I(x)  = I e*^™'  . (67) 

o 

-17  2 

Now,  from  Eq.  (43),  a = 2.0101  x 10  meters  . Since  the  droplet 
velocity  is  assumed  constant,  Eq.  (44)  becomes  ’ 

N 

Av 

where,  from  Eq.  (53), 

N _ 

A q * 

Here,  q is  determined  from  Eq.  (54),  so  that 

“18 

q = 5.424  X 10  coulombs. 
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Then,  since,  .1^  = 96.67  amperes/meters  , therefore, 

15  3 

n “ 1.1052  X 10  particles/meter  . 

Substituting  the  above  values  of  a and  n in  Eq.  (67),  the  approximate 
ion  current  becomes 


I(x) 


^ -0.0222X 

I e 
o 


(68) 


where,  = 8.727  x 10  ^ amperes.  The  approximate  ion  collision  rate 
is  then  obtained  by  substituting  Eq.  (68)  into  Eq.  (41). 


The  dashed  curve  in  Figure  14  shows  the  approximate  ion  collision  rate 
as  a function  of  x,  where  the  collisions  occur  all  along  the  axial  dis- 
tance from  X = 0 to  X.  Note  that  when  x = a,  the  approximate  total  col- 
lision rate  is  only  43  percent  less  than  the  collision  rate  for 
S = 3,000  coulombs /kg. 

In  view  of  the  assumption  (a)  above,  the  probability  density  function 
for  q/m  can  here  be  taken  as  the  normal  distribution  with  a standard 
deviation  of  S = 0. 


5.5  A NON-NORMAL  PROBABILITY  DENSITY  FUNCTION 


5.5.1  niE  CHI-SQUARE  DISTRIBUTION 


The  ion  collision  rate  is  critically  dependent  on  the  probability  density 
function,  f(q/m),  of  the  charge-to-mass  ratio,  as  evidenced  by  Eqs.  (62) 
and  (41).  In  the  foregoing  treatment,  it  was  assumbed  that  f(q/m)  was 
the  normal  (Gaussian)  distribution,  given  by  Eq.  (50).  It  would  be 
interesting  to  allow  f(q/m)  to  assume  a different  form  and  then  learn 
what  effect  this  new  function  has  on  the  ion  collision  rate.  A 
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I 


I 


probability  density  function  which  resembles  a skewed  pormal  distribution 
is  the  Chi-Square  Distribution,  given  by  , 


f -1  2 

£ , (/  s 0) 

2^  T(f) 


where,  v = positive  integer.  Let  z be  a dimensionless^ variable,  given 
by 


2 

where,  C is  a positive  constant.  By  replacing  x with  z in  the  above 
expression,  it  is  easy  to  show  that  the  new  probability  density  func- 
tion for  q/m  is  given  by 


f(q/m)  = 


2^  T(|) 


where , 


J f(q/m)d(q/m)  = 1 


The  mean  and  standard  deviation  of  f(q/m)  are  given,  respectively,  by 


S = ^ (2u)^ 
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Note  that  q/m  and  S are  not  independent,  as  they  are  for  the  normal 
distribution.  We  shall  let  q/m  = 10,000  coulombs/kg.  Hence, 

^ “ TofoW  *^g/coulomb  . 

The  following  value  is  assigned,  v = 4.  Therefore, 

-4 

C = 4.0  X 10  kg/coulomb  . 

Hence,  S = 7071  coulorabs/kg.  Figure  15  shows  a plot  of  f(q/m)  versus 
q/m,  where,  f(q/m)  is  given  by  Eq.  (69).  The  dashed  curve  in  Figure  15 
is  f(q/m),  given  by  Eq.  (50),  where  q/m  = 10,000  coulombs/kg  and 
S = 3,000  coulombs /kg. 


5.5.2  ION  COLLISION  RATE 


From  Eqs.  (69)  and  (47),  the  probability  density  function  of  the  droplet 
radius,  r,  becomes  ' 

V ^ , 


g(r)  = 


2^  T(|) 


ifeY  1 
Ap^J  ^+1 


exp 


-3/2 

2\  2 1^ 

\P  / 


(70) 


From  Eq.  (60),  we  have 

'2  11  V 11/12 

J r^  g(r)dr  J (q/m)  f(q/m)d(q/m)  , 

rj  \ P / (q/m)2 

where,  f(q/m)  is  now  given  by  Eq.  (69),  and,  (q/m)^^  > (q/m)2  . 

The  values  assigned  to  (q/m)^  and  (q/m) 2 should  be  chosen  «io  that  they 
include  practically  all  possible  values  of  q/m.  Since  z(*  C q/m)  obeys 
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Possible  Probability  Density  Functions  for  q/m 


i 

[ 

1 


I, 

I 


I 

i 

f 

i 


the  Chi-Square  distribution,  therefore,  when  u = 4,  99  percent  of  all 
possible  values  of  z lie  in  the  range 

0.208  ^ 14.88  . 

Hence,  99  percent  of  all  possible  values  of  q/m  lie  in  the  range 
520  ^ q/m  ^ 37,200  coulombs/kg  . 


Hence, 


(q/m)j^  = 37,200  coulombs/kg  , 

(q/m)2  = 520  coulombs/kg  . 

Except  for  the  limits  on  the  integrals  in  Eqs.  (58)  and  (62),  these  two 
equations  remain  the  same,  where,  now  f(q/m)  is  given  by  Eq.  (69),  and 
the  integral  limits  are  520  and  37,200.  Employing  the  trapezoidal  rule 
yields: 

•18 

q = 9.77  X 10  coulombs, 

37,200 

J (q/m)  f(q/m)d(q/m)  = 2.537  x 10  ^ (coulombs/kg) 

520 


Hence,  Eq.  (62)  becomes: 


I(x)  = exp 


-6.20(10"^)  4 


e^  ^’‘^f[z(x)]  - 1.01  F(z^) 


(71) 


Substituting  the  above  in  Eq.  (41)  yields  the  ion  collision  rate.  The 
upper  solid  curve  in  Figure  14  shows  the  collision  rate  versus  the  axial 
distance,  x,  from  the  tip  of  the  emitter. 


i 
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5.6  SUMMARY  OF  DISTRIBUTION  RESUr.TS 


The  r<^ults  of  the  foregoing  analysis  can  be  readily  summarized  in 
Figure  14  and  in  Table  VI.  Essentially  five  important  conclusions  can 
be  drawn  from  the  curves  in  this  figure.  They  hre: 


a.  The  most  conspicuous  characteristics  oj  the  four  curves  in' 
Figure  14  is  that  the  collision  rate  ii creases  with  increasing 
values  of  the  axial  distance,  x,  reaching  a maximum  at  the 
termination  of  the  acceleration  region. 

b.  Increasing  the  standard  deviation,  S,  of  the  normal  distribution 
also  Increases  the  collision  rate,  for  my  value  of  x. 

c.  The  shape  of  the  curves  is  the  same  for  >oth  the  normal  (S  ^ 0) 
and  Chi-Square  distributions. 

d.  The  minimum  possible  lower  limit  of  thef integral  in  Eq.  (58), 
and  the  first  integral  in  Eq.  (62),  is  |q/m)2  = 0.  From 

Eq.  (57),  we  see  that,  for  = 10,000|cou' ombs/kg,  this 
minimum  lower  limit  will  occur  when  S ='3333.33  coulombs/kg. 

The  value  S = 3300  coulombs/kg  has  already  been  employed  in 
the  calculations.  From  Figure  14,  we  see  that  the  collision 
rate  for  S = 3300  coulombs/kg  differs,  at  most,  from  the  col- 
lision rate  for  S = 0,  by  approximately  only  an  order  of 
magnitude. 

e.  The  collision  rate  corresponding  to  the  Chi-Square  distribution, 
which  can  be  thought  of  as  a skewed  normal  distribution,  is  only 
slightly  greater  than  the  collision  rate  corresponding  to  the 
normal  distribution,  for  S = 3300  coulombs/kg,  at  every  value 

of  X. 


TABLE  VI 


Probability 
Density  Function 
f(q/m) 

Normal  (Gaussian) 

Normal  (Gaussian) 

Normal  (Gaussian) 

Chi-Square 


COLLISION  RATES 


Mean  Charge- to 
mass  Ratio 
(coulombs/kg) 

10,000 

10,000 

10,000 

10,000 


Standard 

Deviation 

fcoulombs/kf 


3,000 

3,300 

7,071 


Ion 

Collision  Rate 
f collisions /sec) 

1.909  X 10® 

3.327  X 10® 

6.034  X 10® 

7.907  X 10® 
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SECTION  6 

I 

NEGATIVE  ION  PRODUCTION  ' 

6.1  POSSIBLE  MECHANISMS  FOR  NEGATIVE  IODINE  PRODUCTION  ] 

Iodine  exists  in  charged  colloid  droplets  in  three  forms,  Nal,  I , and  | 

I.  Negative  iodine  ions  can  exist  through  ionic  dissociation  of  Nal 

molecules  in  the  glycerol  solvent.  The  presence  of  atomic  iodine  atoms  j 

in  the  droplets  is  probably  due  to  the  solubility  of  iodine  in  glycerol  ] 

after  iodine  ions  have  lost  an  electron  at  the  emitter  electrode.  A 
collision  involving  breakup  or  evaporation  of  a dxoplet  can  release  I" 
directly  by  liberating  the  l“  already  present,  as  indicated  in  Figure  16. 

Atomic  iodine  can  also  be  released  and,  through  a process  of  low  energy 
attachment,  form  I“,  since  iodine  has  a high  electron  affinify.  Attach- 
ment electrons  may  be  provided  from  ionization  events  within  the  droplet 
at  breakup.  Nevertheless,  this  event  has  a low  probal^ility  of  occurring, 
since  electron  capture  only  occurs  at  low  energy  and  ^^he  electrons  gain  | 

high  energy  very  quickly  in  the  intense  electric  fieldy.  ^ 

Negative  iodine  ions  may  also  be  formed  in  collisions  from  the  dissocia-  I 

tion  of  Nal  to  form  Na"^  and  l“.  This  reaction,  as  shown  in  Figure  16,  ‘ 

is  less  probably  than  the  more  energetically  favorable  dissociation  of  I | 

Nal  into  its  atomic  components  of  Na  and  I.  The  latter  reaction  has  a j | 

greater  probability  of  occurring  but,  in  high  energy  collisions  of  I j 

molecular  ions  with  droplets,  more  than  enough  energy  is  available  to  ' ; 

dissociate  Nal  into  its  ionic  components.  f '■ 

I I 

I i 

‘ 1 

Evidence  that  droplets  break  up  in  collisions  can  be  inferred  from  i 

various  optical  and  mass  spectra  data.  The  presence  of  atomic  sodium 
in  the  colloid  thruster  beam  has  been  measured  by  several  groups, 

I ; 

h'  ' ' 

t-  ■ i 

I 
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IN  DROPLET 


OUTSIDE  OF  DROPLET 


I + e-  r 

\ 

OBSERVED 

Na  OPTICAL  SPECTRA 

Na"^  MASS  SPECTRA 
I (OPTICAL  SPECTRA) 


Figure  16.  Mechanisms  of  I”  Production 


Ul 


Including  TRW,  Xerox/EOS,  and  researchers  at  the  University  of  Illinois.' 


Some  evidence  that  atomic  iodine  is  present  in  the  colloid  exhaust  beam 

18 

was  measured  by  Weinstein.  The  latter  study  showed  that  two  and  pos- 


sibly three  iodine  lines  were  obseip^ed  in  optical  spectra  measurements. 
The  identification  of  atomic  sodium  and  iodine  in  the  colloid  beam  poses 
the  question  concerning  their  origin.  It  is  presently  believed  that 
these  atomic  species  are  liberated  from  droplets  which  break  up  in  col- 
lisions with  molecular  ions  or  other  droplets.  Sodium  ions  have  also 
been  detected  using  mass  spectra  technique. 


In  summary,  Na,  I,  and  Na"’’  have  been  observed  in  the  colloid  beam.  It 
is  plausible  to  assume  that  1“  is  also  present  in  the  exhaust  beam.  How- 
ever, l”  has  not  been  measured  and  this  is  partly  due  to  the  difficulty 
of  observing  negative  ion  species  experimentally,  and  before  this  study 
there  was  little  or  no  interest  in  even  considering  i aarching  for  l“. 


6.2  NEGATIVE  ION  TURNAROUND  IN  ACCELERATION  REGION 


Atomic  or  fragment  negative  ions  released  during  a droplet  breakup  event 
will  have  an  initial  velocity  equal  to  the  velocity  of  the  droplet  at 
the  time  of  breakup.  An  illustration  of  the  collision  and  breakup  is 
shown  in  Figure  17.  Whether  the  negative  ion  has  sufficient  energy  to 
escape  the  electric  field  in  the  accelerating  region  will  depend  on  the 
charge-to-mass  ratio  of  the  negative  ion  and  its  velocity  at  the  time  of 


formation. 


Consider  the  one -dimensional  case  of  a charged  particle  in  an  electric 
field.  The  equation  of  motion  for  the  particle  becomes 


dVjxl  _ 0 

a-u-  » 


where  = q/m  is  the  fragment  or  atomic  negative  ion  charge-to-mass 


ratio.  Equation  (63)  determines  the  trajectory  of  the  negative  particle 
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In  the  electric  field  after  its  formation.  It  is  convenient  to  solve 

this  equation  for  the  velocity  in  terms  of  the  potential  at  droplet 

bteakup,  V_,  and  the  potential  at  turnaround,  V , subject  to  the  fol- 
B i 

lowing  boundary  conditions: 


v(x)  = Vg*.  V(x)  = Vg 
v(x)  = 0;  V(x)  = . 


(64) 


The  initial  downstream  velocity  of  the  negative  particle  at  breakup 
becomes 


v/  - 2Cj  <V^  - Vj)  . (65) 

The  velocity  at  breakup,  v„,  can  also  be  determined  from  the  original 
charge-to-mass  ratio  of  the  droplet  which  has  been  accelerated  through 
a potential  difference  Av  = V - V„,  where  V is  the  potential  at  the 

^ O o O 

emitter  surface.  Thus, 

Vg^  = - Vg)  = 2C^AV  , (66) 


where  is  the  charge-to-mass  ratio  of  the  droplet  before  breakup.  If 

we  eliminate  v„  between  Eqs.  (67)  and  (66)  and  solve  for  V , 

B i 


(67) 


Values  for  the  turnaround  potential  were  computed  for  tLo  case  of  the 
pro-'.ictlon  of  negative  Iodine  ions,  assuming  an  initial  droplet  of 

4 

10  r/kg  aixJ  an  »niltter  potential  v^  » 13  kV,  For  this  case. 


- 1 ■» 

1"  arvl  values  for 


are  shown  below  at  corresponding 


VgCkV) 

AV(kV) 

> 

> 

12 

1 

11.98 

11 

2 

10.97 

8 

5 

7.93 

5 

8 

4.89 

3 

10 

2.86 

1 

12 

0.84 

These  results  Indicate  that  negative  iodine  particles  do  not  continue 
very  far  after  their  formation  in  the  accelerating  region.  At  the  boun- 
dary of  the  accelerating  region  where  V(x)  = = 0,  we  can  solve  Kq.  (67) 

for  AV  to  obtain 


For  iodine  negative  ions  we  find  AV  = 12.83  kV.  Therefore,  iodine  ions 

formed  after  a droplet  lias  been  accelerated  through  a potential  differ- 
ence greater  than  12.83  kV  will  escape  tlie  accelerating  region  (V^  = 13  kV). 
We  conclude  from  this  analysis  that  nearly  all  the  negative  iodine  ions 
formed  within  the  accelerating  region  do  notescape  and  will  be  accelerated 
toward  the  emitter. 


I 
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SECTION  7 

CONCLUSIONS  AND  RECOMMENDATIONS 


7.1  SUMMARY  OF  THE  MODEL 

This  report  has  presented  the  essential  features  of  a model  describing 
the  mechanisms  of  colloid  emitter  damage  responsible  for  the  degradation 
of  performance.  The  model  is  summarized  here  and  the  reference  section 
listed. 

Section 


a. 

Particle  species  in  the  exhause  charged  droplets 

(3) 

Charged  droplets 

Molecular  ions 

Neutral  molecules 

Electrons 

b. 

Highest  density  of  particles  in  the  portion  of 
the  exhaust  beam  within  the  acceleration  region 

(3) 

Molecular  ion  current  density 

S3  6 X 10^^  ions/sec 

(3.2) 

9 3 

Charged  droplet  density  S3  10  /cm 

(3.3) 

; 

9 3 

Neutral  molecules  =al0  molecules/cm 

(3.1) 

Electrons,  very  few  expected. 

c. 

Charged  droplet  distribution  is  represented  by 
the  mean  charge-to-mass  ratio. 

(3.4) 

i ; 

d. 

The  mean  droplet  has  a total  molecular  binding 
energy  of  3 keV. 

(3.4) 

1 

I' 

e. 

The  acceleration  region  also  has  a highly 
divergent  and  intense  electric  field. 

(2.2) 

1 . 

i 

f. 

Collisions  will  occur  between  these  particle 
species  because  they  travel  with  greatly 
different  velocities. 

(4) 

1 
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p 

1 

t 

1 

i 

1 

i 

1 

Section  j 

g> 

t 

Only  collisions  between  charged  droplets  and 
both  molecular  ions  and  electrons  could  cause 
droplet  breakup. 

(4) 

h. 

Molecular  ion  collisions  with  charged  droplets 
occur  at  the  higliest  rate  because  of  the  current 
densities  and  collision  cross  sections.  These 
rates  range  between  10^  and  10^/sec.  Inclusion 
of  small  droplet  collisions  with  large  droplets 
increases  the  above  rate. 

(4.3) 

1. 

Droplet  fragment  particles  have  been  observed 
experimentally,  providing  evidence  of  collisions. 

(4.1) 

j. 

Charged  droplets  may  decrease  in  size  due  to 
evaporation  and  result  in  positive  charge  emission. 
This  is  not  likely  to  lead  to  negative  ion 
production. 

(4.4) 

k. 

Several  distributions  of  particle  sizes  were 
introduced  into  the  model  and  compared  with  the 
single  particle  model. 

(5) 

1. 

• 

Normal  distributions  having  various  standard 
deviations  and  also  a Chi-square  distribution 
were  used  with  no  significant  alteration  of 
the  results. 

1 

(5) 

m. 

When  molecular  ions  collide  with  droplets,  the 
droplets  break  up  due  to  the  deposition  of  the 
ion  energy  in  the  droplet. 

(4.3) 

n. 

Several  processes  are  advanced  for  the  produc- 
tion of  negative  iodine  at  the  collisional 
breakup  of  a droplet. 

(6.1) 

o. 

All  negative  ions  (I")  produced  in  the  accel- 
eration region  from  the  breakup  of  droplets, 
except  for  a few  near  the  edge  of  the  accel- 
eration region,  do  not  escape  from  the  region. 

(6.2) 

P* 

These  negative  ions  are  returned  to  the  emitter 
and  bombard  the  surface  with  energies  of  several 
thousand  electron  volts. 

( 

(6.2) 

q* 

Each  ion  sputters  between  1 and  10  atoms  from 
the  emitter  surface,  depending  upon  the  bom- 
barding ion  energy. 

( 

(2.3) 

r. 

* 

Tlie  observed  erosion  is  adequately  explained  by 
this  model  if,  on  the  average,  each  molecular  ion/ 
droplet  collision  results  in  the  production  of  one 
negative  ion. 

f 

(2.3)  , ! 

. J 
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Section 


s.  The  eroded  emitter  surface  causes  a change  in 
the  intense  electric  field  at  the  emitter. 


t. 


Alteration  of  the  field  at  the  emitter  results 
in  degradation  of  performance. 


(2.1) 

(2.1) 


7.2  SUGGESTED  EXPERIMENTS  FOR  DETECTION  OF  NEGATIVE  IONS 


Investigations  into  the  possible  mechanisms  of  negative  ion  production 
suggest  several  experimental  directions.  Understanding  of  the  colloid 
emitter  erosion  process  may  require  a two-phase  experimental  approach. 
One  experimental  phase  would  treat  the  problem  of  detecting  negative 
ions  and  identifying  the  species.  Measurements  using  mass  spectrometer 
techniques  will  yield  this  type  of  information.  Other  experiments 
could  be  designed  to  provide  information  concerning  the  mechanism  of 
negative  ion  formation. 


Two  experimental  directions  have  been  considered  to  provide  data  that 
could  verify  models  of  negative  ion  formation.  One  such  experiment 
woi^d  involve  a cross-beam  technique  using  two  colloid  sources.  One 
source  would  be  operated  in  the  molecular  ion  emission  mode  and  one 
operated  in  the  heavy  charged  droplet  mode.  If  I”  or  other  negative 
ions  are  released  during  droplet  breakup  collisions,  absorption  spectra 

?y  provide  a suitable  detection  method.  A less  direct  approach  for 
sting  the  hypothesis  of  negative  iodine  liberation  in  droplet  colli- 
sions involves  absorption  spectra  measurements  in  the  region  of  a 
f^ollector  downstream  from  a colloid  source.  Droplet  impact  at  the 
collector  involves  higher  collision  energies  than  are  generally  avail- 
able within  the  accelerating  region,  but,  nevertheless,  may  provide 
useful  clues  as  to  the  formation  and  type  of  negative  ions  produced 
when  charged  droplets  break  up.  In  both  droplet  breakup  experiments 
suggested  above,  mass  spectrometer  measurements  could  be  substituted 
for  absorption  spectroscopy  to  detect  negative  ions. 
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7-3  R!  :OMMENnATIONS 


Tc'  reti'.;  ; emitter  erosion  and  the  subsequent  performance  degradation, 
seVerc.  i ecommendations  were  made: 


I 


i 


- Operation  at  lower  voltages  would  reduce  the  bombardment 

energies  of  the  negative  ions  and,  therefore,  the  sputtering 
rate. 

b.  Use  a propellant  dopant  having  a much  lower  mass  than  iodine 

to  reduce  the  sputtering  rate.  i ; 

, t 

t Use  an  emitter  material  having  a lower  sputtering  yield.  ( * 

Use  a propellant  containing  no  atoms  or  molecules  having  a 
significant  electron  affinity,  thus  reducing  the  production 
of  negative  ions. 

1 Eliminate  the  productivity  of  positive  ions  at  the  emitter 
which  therefore  reduces  the  ion/droplet  collision  rate. 

Elimination  of  positive  ion  production  also  increases  the 
efficiency  of  the  thruster  operation. 


I 

i 


I 

I 
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